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OPTICS OF METALS 


The following four contributions were read at the Society’s meeting at 
Southampton on 18th December 1951. 


A New Photographic Method of Measuring the Dispersion 
of the Optical Constants of Metals 


By J. BOR 
West Ham College of Technology 


MS. received 2nd April 1952; read before the Society at Southampton on 18th December 1951 


ABSTRACT. ‘The information from which the optical constants are determined is 
produced on a photographic plate as a vertical band of variable intensity for each wavelength. 
‘This is compared microphotometrically with an adjacent band of graded intensity produced 
by the same wavelength. 


measurement of the ellipse of vibration produced when plane polarized 
light at an azimuth of 45° is reflected obliquely from a plane metallic 
surface. ‘This requires the simultaneous determination of two independent 


T= basic principle of this method, as in most polarimetric methods, is the 


Fig, 1, 


quantities, and for this reason it has hitherto been possible to measure only one 

wavelength at a time. The most usual quantities are the ratio of the vertical 

and horizontal components ¢% (fig. 1), together with the phase angle A between 
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them or, as used here, the ratio of the major and minor axes of the ellipse y and 
the orientation of the major axis to the vertical 7. The method of calculation 


used is first to transform from 7 and y tos and A by the geometrical transformations _ 


cos 2y;=cos 27 cos 2y, tan A= —tan 2n cosec 2y, and then to calculate the 
optical constants from the equations 
2 sin? ¢ tan? d (cos A+sin 2) 
sin 2%(cos A+cosec 2)?’ 
sin? ¢ tan? ¢ sin A 
tan 2x%(cos A+cosec 2x5)’ 
where ¢ is the angle of incidence, m the refractive index, k the extinction 


coefficient, « the dielectric constant, o the electrical conductivity at frequency », 
7 the period of the light, i.e. t=1/v. 


«=n? —k* =tan’ ¢ — 


ol=nk = 


The dispersion curve of the optical properties of metals is generally accepted | 


as being of much greater importance than very exact measurements of one 


. . . { 
wavelength only. ‘To make measurements over a wide spectral range it is 


| 


necessary to observe one wavelength at a time over a lengthy period during which 


the optical properties, particularly of films, may change considerably. The 
apparatus to be described (Bor and Chapman 1949) enables all the information 
necessary for the calculation of the optical constants between approximately 
0:41 to 0-651 to be recorded on a single photographic plate in the form of a 
pattern of variable density. 


Cylindrical 


7 RE: Gea) [E 


5 Newel cae 
Rotator Analyser L, tall 
Camera 


Kige2s 


The optical system is shown in fig. 2. Parallel reflected elliptically polarized 
light is passed through the rotator, the analyser, a dispersing system of two prisms 
(one only shown in the diagram) and then collected by the two cylindrical lenses 
L,, L, to form a pattern on the photographic plate P. The rotator consists of 
two prisms of right- and left-handed quartz in optical contact. Above the céntral 
plane the ellipse will be rotated clockwise, below counterclockwise. Figure 3 (a) 
indicates the relative orientations of the ellipse in the light emerging from the 
rotator. 

If CC is the plane of no rotation, the orientation of the ellipse is given by 
7 =(AC/AB)z. This central plane is defined by a very fine wire W (fig. 2) placed 
immediately after the rotator. A magnified image of the wire and of the pattern 
of light emerging from the rotator modified by the analyser with its vibration 


direction vertical is thrown by the cylindrical lens L,, of focal length 20 cm at — 
its conjugate focus, on to a photographic plate P, the refraction taking place in a 
vertical direction oni Refraction in the horizontal direction only is produced | 


by the lens L, of focal length 65 cm (the length of the camera), and this converges 


the parallel light i in the horizontal plane on to the photographic plate P. This 
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combination of lenses produces a narrow vertical band of intensity which is not 
uniform (fig. 3 (b)), but which oscillates in the ratio a2/b2. Hence the square root 
of the ratio of minimum to maximum intensity gives the ratio y of the axes of 
the ellipse. Dispersion of the white light by a train of two prisms produces a 
spectrum crossed by a system of bright and dark bands together with the fiducial 


After 
Analysis 


(a) (0) 
Fig. 3. 


After Rotation Intensity 


line. Adjustment is provided so that the photographic plate can be suitably 
tilted to allow for variations in focal length of the different wavelengths. The 
cylindrical lenses can be moved horizontally and can be independently rotated 
by slow motion about a horizontal axis. Accurate adjustment is obtained by 
using the line spectrum from a mercury arc, lens L, being adjusted before L, is 
placed in position. The adjustment about the axis is very sharp, a very small 
rotation of either lens producing an indistinct image. 


|UD UTA 


My 


TOIT 


Fig. 4. Double triangular slit. 


In order to translate relative blackenings on the photographic plate into true 
relative light intensities the following procedure was adopted. ‘The rotator and 
analyser, which are mounted together as one unit, are removed and replaced by an 
accurate double triangular aperture (fig. 4) arranged vertically, so that its plane 
passes through the position previously occupied by the fiducial line. This 
produces a spectrum in which, owing to the action of the cylindrical lenses, the 
light intensity for each wavelength increases linearly with distance measured 
vertically from the centre. ‘The two spectra are interleaved by placing a grid G 
(fig. 2), of alternate equal opaque and transparent strips, in front of the photographic 

3 D-2 
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plate. An exposure of about 5 seconds is first made with the rotator unit in 
position. The latter is then removed and replaced by the double triangular slit, 
the photographic plate is moved horizontally the appropriate distance in its own 
plane, so as to uncover the unexposed segments of the plate, and another exposure 
of about 5 seconds is made. ‘Thus adjacent strips of the two spectra are formed 
by the same wavelength. Figure 5 shows the pattern obtained with aluminium, 
the whole operation taking about 2 minutes, as compared with the several hours 
necessary in covering the same spectral range one wavelength at a time. - 


0-4 ps ; : 0:65 uw 
Fig. 5. Light pattern obtained with aluminium. 


0-4 w 0-65 uw 
Fig. 7. Calibration light pattern. 


Light and dark are reversed on the photographic plates. 


A self-luminous point source of light collimated by a lens or mirror would 
be the ideal illuminant, since it would produce a complete pattern, and would 
be of uniform intensity over the field of the rotator. Since in practice a light 
source of finite size must be used, its permissible dimensions are governed by the 
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following considerations: (a) The width must be such that the wavelength 
confusion in the spectrum due-to horizontal divergence of the beam must be 
sufficiently small compared with the variation of optical constants with 
wavelength. (6) The permissible length is governed by the maximum vertical 
divergence of the beam. This produces confusion of the fringes because of the 
double refraction of a light ray which passes through quartz at a small angle to 
the optic axis. Calculation showed that a divergence of about ++4° would 
produce a negligible confusion. 

These conditions were fulfilled by using a headlamp bulb with a small, 
straight, compact filament of size about 4 x 1 mm at the focus of a large aperture 
50 cm achromatic lens. ‘The wavelength confusion was no more than 25A in 
the blue and 1004 in the red, and the uniformity as tested by a photocell along 
a spectrum line (rotator unit and aperture removed) was quite constant and 
varied by only $°% at the extreme ends, which were avoided when measurements 


‘were made. 


A Hilger non-recording microphotometer fitted with a 25 cm horizontal 
scale reading by vernier to 0-01 mm was used for measuring the intensity of the 
light transmitted through the plate and for fixing the positions of the maxima 
and minima. 

Measurement of ratio of axes. ‘Transmitted light intensities at the positions. 
of maximum and minimum are matched with the adjacent graded intensity. 
Since maximum incident light on the photographic plate produces maximum 
blackening and minimum light transmission it follows that the ratio of the major 
to the minor axis is given by (L,/L,)"? (fig. 6). Mean values are taken for all 


Double Slit : 
Intensity Fringe 
Intensity 


Transmitted Intensity — 


r Lp “4 
Fig. 6. Transmitted intensity distribution through the photographic plate. 


the maxima and minima for each wavelength. It is preferable to keep within the 
linear portion of the characteristic curve of the plates, which were Ilford special 
rapid panchromatic developed to a y value of about 1. ‘The linear region extends 
over a density range of 0-4 to 2:0, covering an intensity ratio of about 60:1, 
equivalent to an axial ratio of 7-7 to 1. For the metals measured the axial ratio 
did not exceed 3:1. 

Measurement of orientation. ‘The position of the fiducial line relative to 
the central plane was measured on the calibration plate (fig. 7) by allowing plane 
polarized light vibrating in a horizontal plane to pass through the rotator and 
analyser, adjusting the wire to be parallel to the central minimum of transmission 
which coincides with the position of no rotation, and then measuring its distance 
from this line. The distance between the fiducial line, marked by absence of 
blackening, and the position of the various maxima and minima on the specimen 
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plate (fig. 5) are measured, from which mean values of AC and AB (fig. 3(a)) 
are obtained. In practice, readings of equal deflection are taken on both sides 
of the various maxima and minima. 

The apparatus was tested by measuring a specimen with the author’s 
photoelectric apparatus (Bor 1937, Bor, Hobson and Wood 1939) and then with 
the new apparatus. The results agreed to within the limit of experimental error 
throughout the spectrum range. 

Measurements have been made im vacuo on gold, copper and aluminium 
films prepared by evaportiens and the results will be reported in a subsequent 
communication. 

REFERENCES 
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Photoelectric Analysis of Elliptically Polarized Light 
By J. F, ARCHARD, PP -CLEGG AND Ay Oia ior 


University of Southampton 
MS. received 21st March 1952; read before the Society at Southampton on 18th December 1951 


ABSTRACT. Two photoelectric methods of analysis of elliptically polarized light have 
been devised. One is a development of that of Kent and Lawson in which a rotating 
analyser yields a zero signal only when the light is circularly polarized. We use a double 
image prism after a rotating analyser and the two resulting modulated beams are in anti- 
phase only when the axes of the double image prism coincide with the axes of the vibration 
ellipse. In our second method, the -light is modulated by chopping, and a double image 
prism selects two components a phase. In principle these are equal only when the axes of 
the double image prism are at 45° to the axes of the vibration ellipse. Thus the directions 
of the major and minor axes are found, and if the outputs from the two photocells which 
receive the two components are balanced on an electrical bridge the axial ratio may be 
obtained. The double image prism plays here the part of the half shadow device in visual 
observations, but as its half shadow angle is 90° instead of only a few degrees, the photocells 
are used under optimum conditions, i.e. the signal-to-noise ratio is maintained as large as 
possible and the potential accuracy of the method greatly exceeds the limitations imposed by 
the rest of the apparatus. By slight modifications the method may be made partially 
self-recording. 


Some examples are given of measurements made on metallic films undergoing surface 
changes. 


§1. INTRODUCTION 

RANSMISSION through birefringent material, or reflection at an interface 

| causes linearly polarized light to become elliptically polarized. Measure- 
ment of the effect is not identical in the two cases. When birefringence 

is not accompanied by perceptible dichroism, the phase difference A between 
the two components of light parallel and perpendicular to the two relevant axes 
of birefringence is the only subject of measurement. If dichroism is present, 
these two components suffer differential absorption, and it is necessary also to 


find the ratio of these components before and after transmission by the material. 
If, relative to the axes of birefringence, y be the azimuth of the incident plane 


lari : 
polarized light, then T tan x =tan x’ eeeeee (1) 
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where T is the ratio of the transmission coefficients, and tan y’ is the measured 
ratio after transmission. 

In the case of reflection experiments, the reference planes are normal to the 
incident or reflected wave front, and parallel and perpendicular to the plane of 
incidence. ‘The problem is then similar to the simultaneous determination of 
birefringence and dichroism, the two quantities of interest now being A, the 
differential phase change between the two components introduced on reflection, 
and p the ratio of their amplitude reflection coefficients, where 


tan y’ 
Re § 
an x 


=tan ¢. RTs 


us is the so-called azimuth of restored polarization. 

In what follows, we shall be concerned primarily with the reflection case, 
but the methods described could equally well be applied to the transmission 
problem. 

1.1. Visual Methods 


When birefringence alone is in question, measurement of the phase 
difference A is usually made by means of an optical compensator. ‘The principle 
of some, such as that of the Babinet or Soleil type, is that the compensator 
introduces a phase change equal and opposite to that introduced by the 
specimen; the elliptically polarized light is thus reduced to plane polarized light 
with the same azimuth as the original incident light. 

A compensator much used is that known as the Senarmont, which consists 
simply of a quarter-wave (\/4) plate. One method of using it, which has 
similarities to the photoelectric methods discussed later, consists in allowing y 
and consequently y’ to have an arbitrary value, and adjusting both the Senarmont 
compensator and the analyser until extinction is obtained. When this is done 
the axes of the compensator coincide with the major and minor axes of the 
vibration ellipse, from which the direction of these axes is known, and from the 
position of the analyser the ratio of the axes is determined. If y be the azimuth 
of the compensator, and B+y the azimuth of the emergent plane polarized 
light, it may be shown that 


cos 2’ =cos 2y cos 28, 
tan A=tan 28 cosec 2y. 


Thus A the phase change and y’ can be calculated. The ratios of the 
transmission or reflection coefficients are determined by eqns. (1) or (2), from 
the calculated value of y’, and from the azimuth of the incident plane polarized 
light y. 

In the above arrangement there are three variables y, y and 8. It is possible, 
and sometimes convenient, to fix one of these three. The procedure sometimes 
adopted is to set y at +45°, i.e. in the case of reflection the azimuth of the 
incident plane polarized light must be set so that the axes of the reflected 
elliptically polarized light are at 45° to the plane of incidence. ‘This arrangement 
makes it possible to determine 4 and A directly, since in this case 


eee oe (4) 


Physically, these equations mean that the polarizer is set so as to compensate 
for the difference in the reflection coefficients, the two components (parallel and 
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normal to the plane of incidence) being equal after reflection. The axes of the 
reflected ellipse are then at 45° to the plane of incidence, and their amplitude 
ratio is tan $A. 


1.2. Photoelectric Methods 


Analysis of elliptically polarized light has been made in non-visible regions. 
of the spectrum by ingenious photographic methods such as those due to Voigt 
(1901) and Minor (1903), to Szivessy, Dierkesmann and Munster (1933) and to” 
Bor and Chapman (1949), but in no case is the method very rapid, for the 
exposure may take up to very many minutes, and much work has to be done in 
measuring the parameters of the patterns obtained. Sometimes this is not an 
insuperable objection, but if it is wished to follow the change of ellipticity as 


the surface undergoes physical or chemical change, then a faster method must _ 


be devised. The use of photoelectric or electronic devices seems obvious, but it 
appeared advantageous not to adapt the conventional visual methods to suit 
photoelectric detection, but to devise methods new in themselves. 


Any photoelectric method should use to the best advantage the properties of — 


the photoelectric cell. At low intensities of illumination it has little advantage 
over the eye, but for higher intensities it is greatly superior. A number of 
_workers, Bruhat and Guinier (1933), Rank, Light and Yoder (1950), and Levy, 
Schwed and Fergus (1950), have applied photoelectric cells to polarimetry, and 
Bruhat and Grivet (1934) to the analysis of elliptically polarized light. Broadly 
speaking, these attempts, though ingenious, have not made the best use of the 
properties of the photocell as they have generally merely replaced the eye by a 
photocell in existing optical apparatus. 

A promising method is that described by Kent and Lawson (1937), who. 
used the fact that rotation of an analyser produces no change in intensity when 
the light is circularly polarized. This method appears very suited to the 
characteristics of a photocell, e.g. ‘dark current’ is of no importance, since the 
rotating analyser modulates the light beam, and the photocell is used as an 
a.c. device. Moreover, a narrow band amplifier used between the receiver and 
detector gives a greatly improved signal-to-noise ratio. We have examined 
this method in more detail and, as discussed below, have found a number of 
disadvantages which are not immediately obvious from the account given by 
Kent and Lawson. Moreover, the method as described by these authors must 
employ some form of optical compensator if it is desired to detect elliptically 
polarized light; the elliptically polarized light must be transformed by the 
compensator, not into plane polarized light as in visual methods, but into 
circularly polarized light. 

We have, therefore, considered the possibility of analysing elliptically 
polarized light without the use of an optical compensator, the orientation y and 
ratio tan f of the axes of the ellipse (see eqns. (3) above) being determined directly. 

‘T'wo such methods, previously described in a preliminary report (Archard, 
Clegg and ‘Taylor 1950), are discussed in more detail below. The essential 
technique involved is that a double image prism is used to select from the 
elliptically polarized light two components polarized in mutually perpendicular 
directions; each of these components is received by a separate photocell and, 
as the light is modulated, either mechanically or optically, the resultant signals 
may be matched with one another. Depending on the method of modulation, 
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these signals may be in phase (Method A), or in anti-phase (Method B), but in 
both cases the matching circuit is essentially the same. 
The type of matching circuit used is shown in fig. 1. The two signals are 


_ fed through cathode followers V, and V3, acting partly as buffer stages, to 


opposite halves of the twin triode matching valve V,. From this stage the output, 
zero when matched, is fed to a cathode-ray oscillograph via a tuned amplifier 
(Sturtevant 1947). Helical potentiometers R, and R,, giving 3600 degrees of 


‘* [ * 220v 


haduaa C4 


cooccoasccsom 


Fig. 1. 


spindle rotation, in the cathode circuits of the cathode followers, enable the 
signals to be matched and their ratio determined. ‘The output of the matching 
circuit is fed to the tuned amplifier from either A or B, depending on whether 
it is desired to balance in- or out-of-phase signals. 


1.3. General Discussion of Sensitivity 


The sensitivity of any photoelectric method will depend primarily on the 
signal-to-noise ratio. Noise arises chiefly as diode noise in the photocells, and 
Johnson noise in the grid resistors of the first valves. We shall make the 
following general assumptions: intensity of elliptically polarized light 
J, =10~* lumen (mercury green light was used in our experiments); photocell 
sensitivity =30 wa/lumen; grid resistors R=50 megohms; amplifier band- 
width = 10 c/s. 

On this basis, the noise at the grids of the first valves consists of about 4 uv 
of diode noise, and 3 pv of Johnson noise, or a total of about 5v. If it is assumed 
that the minimum signal-to-noise ratio should be about 4, then it follows that 
the minimum detectable signal is equal to that produced by a variation of light 
intensity J =1-3 x 10-* lumen, i.e. 62/Jy)=1-3 x 10-4. 

This criterion will be used to derive the sensitivity of the methods discussed 
below. It may not, in fact, be universally applicable, but it will suffice to show 
that under most conditions these methods are extremely sensitive, and compare 
very favourably with the usual visual methods. 


§2. THE ROTATING ANALYSER 
We have applied the rotating analyser method of Kent and Lawson to the 
detection of circularly polarized light reflected from the metallic surface, the 
incident light being plane polarized. The arrangement used is shown schematically 
in fig. 2. The analyser was rotated at about 43 rev/sec, and the amplifier tuned to 
86c/s. The angle of incidence @ and the polarizer azimuth were adjusted until 
no detectable output was observed on the cathode-ray oscillograph. ‘This occurs 


at the principal angle of incidence @ (the angle of incidence at which the differential 
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phase change A is 7/2). There are four positions of the polarizer for zero output, 
two pairs each separated by an angle 2%, where y is the principal azimuth (the 
special value of y;when 6=6). Thus, applied to a solid metal surface, the method 


enables 0 and y to be determined directly, and from these angles the optical 
constants of the metal can be determined (O’ Bryan 1936). 


Tuned 
Amplifier 


4 Plate 
Photocell 


Rotatin 
Analyser 


Fig. 2. 


On the signal-to-noise considerations given above, the theoretical sensitivity 
of the adjustments corresponds to errors in # and A of about ten seconds of are. 
In fact, the instrumental accuracy of the apparatus was only of the order of one or 
two minutes and, as discussed below, other considerations affect the final accuracy. 
However, the fact that internally consistent results can be obtained may be judged 
from four sets of consecutive readings made on an aluminium surface (table 1). 


Table 1 
Run no. 1 2 3 4 
6 75° 287 75° 30’ 75° 29’ 15° 20 
yb 40° 50’ _ 40° 54’ 40° 55’ 40° 52’ 


2.1. Sources of Errors 

Consider the effect of small misadjustments from their null positions of the 
polarizer azimuth or of the angle of incidence (this latter is effectively a change in A) 
for the arrangement shown in fig. 1. When both adjustments are in their null 
positions circularly polarized light is obtained, and ideally there is no output signal. 
If one of the two adjustments is moved slightly, an ellipse of small ellipticity results, 
and a small output signal is obtained. ‘The orientation of the ellipse and the phase 
of the signal differ according as to which of the two adjustments is altered. The 
effects may be summarized as follows : ; 


Misadjustment Orientation of ellipse Output Signal 

Polarizer 45° to plane of incidence E=+E, cos2wt 

Incidence parallel or normal to plane of E= + E,sin2ot 
incidence 


where 27rw is the frequency of rotation of the analyser. 

The sign of the output signal is dependent on the sign of the error of adjustment. 
Thus polarizer settings of +(%+6x) and —(%+8x) produce identical signals. 
Alternatively, for a correct setting of the polarizer azimuth but a small misadjust- 
ment of the incidence, the phase of the signal depends on which of the four 
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positions of the polarizer azimuth is adopted. Thus a given misadjustment of 
incidence produces a signal + FE sin 2wt at one polarizer position, and — E sin 2wt 
at the other. 

‘These considerations have an important bearing on the detection of errors. 
A number of causes may produce variations in the light flux, or signals from the 
photocells, which have a component at the required frequency (86 c/s in our Case). 
The adjustments to the apparatus now have to be made so that the required signals, 
due to the polarimetric effects being studied, are equal in amplitude and opposite 
in phase to the unwanted signals. The implication of the last paragraph is that the 
cos 2wt term of such unwanted signals will produce errors in the polarizer readings 
which are not eliminated in a complete set of readings. On the other hand, the 


sin 2wt component will produce errors in of +80 and — 80, for polarizer azimuths 


of +x and —Y respectively. These errors in #, which have been observed in 
practice, are averaged out in a complete set of readings. 

The following causes have been shown to produce errors of this type : (@) 
vibration, (6) transmission variations in the rotating analyser, (c) stray light 
reflected from the rotating analyser assembly, (d) the selective photoelectric effect 
at the photocell, (e) angular deviations of the light beam produced by the rotating 
analyser. 

It has proved extremely difficult to obtain calcite polarizing prisms which do 
not produce a detectable deviation of the light beam. Apart from the special 
consideration of unwanted signals, this deviation automatically reduces the 
accuracy of the measurements, because the diameter of the exit pinhole must be 
increased. 


Chopper 


ae Polarizer aL ea 3 


Auxiliary 
Analyser 
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in-phase Amplifier 
Signals 


Fig. 3. 


The selective photoelectric effect, combined with the semi-cylindrical shape of 
-the photocell cathode, means that the sensitivity of the photocell is dependent 
-on the azimuth of incident plane polarized light. In our case, the effect was 
-such that quite small movements of the photocell made it necessary to change the 
polarizer azimuth by as much as one degree. ‘The effect was eliminated by fixing 


_aA/4 plate in the same mounting as the rotating analyser, with its axes set at 45° to 


the plane of polarization of the analyser. This converted the plane polarized light 


-emergent from the analyser into circularly polarized light. 


§3, DIRECT ANALYSIS OF ELLIPTICALLY POLARIZED LIGHT: 
METHOD A 
The method, an experimental arrangement of which is shown schematically in 
‘fig. 3, uses a double image prism and is similar to a visual method previously 
described by Ingersoll and Littleton (1910) and a photographic method due to 
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Pfestorf (1926). The double image prism selects two mutually perpendicular 
components of the elliptically polarized light, and the principle involved may be 
simplified as follows : (a) When the two components are equal, the double image 
prism is set at 45° to the axes of the ellipse : this determines the orientation of the 
ellipse. (b) The double image prism is then set along the axes as determined in 
(a), and the ratio of the components gives the intensity ratio of the axes of the 
ellipse. ‘The ellipse is thus completely determined, except for its sense of 
rotation. 

The method of using the apparatus will now be described in more detail. In 
general, no matter what the setting of the scaling potentiometers, there will be 
four positions of null output as the mounting is turned through 360°. ‘These 
positions will be in two complementary pairs lying symmetrically with respect to 
the axes of the vibration ellipse. By suitably adjusting the bridge potentiometers, 
these null positions can be obtained at about 45° to the ellipse axes, and this can 
be shown to be the most sensitive position. By taking the mean value of these 
four positions, the direction of the axes of the vibration ellipse may be 
determined. 

Once the axes have been located, the intensity ratio along them may be found by 
turning the mounting until the axes of the prism lie along those of the ellipse. Ifthe 
measured ratio is p, in this position, then when the mounting is turned through 90°, 


the new ratio will be p,, the change being due to the differences between the two ~ 


circuits (amplifications, photocell sensitivities, etc.). The true ratio is given by 
(Pi p>)”, and the amplitude ratio, the desired quantity, by (p,p.)!*._ In practice 
the ratios p, and p, are each found by taking the mean values of the ratios measured 
at positions 180° apart. 

The angles y and f are thus determined, the first directly and the second from 
tanB=(p,p.)"*. yx’ and A the required parameters of the ellipse may then be 
calculated from eqns. (2) and (3). 

If the mounting be turned so that the axes of the prism lie in and perpendicular 
to the plane of incidence, then by setting the bridge fora null, tan y’ may be measured 
from this one reading only. ‘This is true, provided that the bridge ratio corres- 
ponding to equal intensity in the two beams from the prism be known. This 
“equality ratio’ may be found, for example, by measuring the ratio for a null with 
plane polarized light incident on the prism in azimuth 45°. 

By analogy with the use of the Senarmont compensator A and % may be 
measured directly if the polarizer azimuth x is adjusted till y =45°, their values being 
calculated from eqn. (4). 

Provided that a Wollaston rather than a Rochon double image prism be used 
(Archard 1950) it would be possible to match the signals by means of an auxiliary 
analyser, as is indicated in fig. 2, instead of by adjusting the bridge as has been 
described. ‘This analyser, an ordinary polarizing prism, should be placed immedi- 
ately after the double image prism, preferably in the same mounting so that its 
orientation with respect to the double image prism can be read directly. When 
the double image prism is aligned along the ellipse axes, the axial ratio can be 
measured from that orientation of the analyser which results in a null output. 

If the electrical circuits associated with the two photocells were adjusted for 
approximate equality, the average value obtained from the four positions of the 
auxiliary analyser would suffice to give the axial ratio. This analyser could be set 
at 45° in the determination of the orientation y by the procedure described above. 
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- 3.1. Sensitivity 


Assume that the method described above has been set up, using an auxiliary 
analyser. It is required to measure the orientation y and ratio tan f of the axes of 
the ellipse, these being given directly from the settings of the double image prism 
and auxiliary analyser respectively. The signals resulting from angular errors of 
adjustment dy and 68 may be regarded as being due to light intensities 6, and 
5I, respectively. Then, it may be shown that 


Bins cos 2B5y, blogs SIZ pop ste Pero OC, (5) 
Ih Io 
where J, is the intensity of the elliptically polarized light being analysed. 

Once again, applying the criterion of § 1.3, it may be seen that undera wide range 
of conditions the theoretical sensitivity in the measurement of y and 8 is much 
less than one minute of arc. 

Usually the quantities required are % and A, and therefore it is necessary to 
convert these errors in y and f to errors in # and A, but, except under the most 
unfavourable conditions, they will be of the same order. In the special case where 
y is set at 45°, it has been shown that A=28, and ¢ is given directly from the 
polarizer setting. In this case, the errors in # and A are given by the equations 


eee Ait ein ASN) $0 Paes (6) 
I Io 
3.2. Results 


Though, as has been shown, the method is very sensitive, a high experimental 
accuracy has not yet been attained, or even attempted. The angular scales in our 
apparatus were divided only to tenths of a degree, which gave sufficient accuracy for 
present purposes. The internal consistency of the results, a measure of the sensi- 
tivity of the method, was good, as can be seen from table 2. The results refer to 
values for light reflected from a thin gold film, plane polarized light at 45° azimuth 
being incident at 45°. 


Table 2 
y 24° 40’ 24° 37’ 24° 38’ 24° 39’ 
B 127 05% 1052 L2S0 54 1272067 
if 26° 46’ 26° 43’ 26° 44’ 26° 45’ 
A 307 36" 30° 39’ SO? Sek AO? Bek 


Much work has been done with the apparatus in determining the values of 
and A as a function of thickness for thin films of various metals evaporated and 
measured in vacuo (Clegg 1952). As an example of measurements within the 
scope of the apparatus, fig. 4 shows the variation of s and A with thickness for a 
copper film; the thickness was determined by means of the microbalance (Clegg 
and Crook 1952), the density of the film being assumed to be that of bulk copper. 
Figure 5 shows an automatic recording of the oxidation of a thin copper film which 
was obtained as follows: after deposition of the copper film the bridge was balanced 
so as to measure tan y’ directly as has already been indicated; air was then admitted 
to the system and the film partly oxidized, thus altering the value of y’. The 
bridge, which was not touched, was then no longer balanced, and it was this 
out-of-balance signal that was suitably recorded. 
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§4. DIRECT ANALYSIS OF ELLIPTICALLY POLARIZED bIGH is: 
METHOD B 


In Method A, described above, the bridge is balanced for signals in phase, and 
two separate adjustments of the apparatus are necessary to obtain the orientation y 
and amplitude ratio tan of the axes of the ellipse. ‘The method described below 
utilizes the fact that the bridge must be balanced both in amplitude and phase. In 
this way it is possible to obtain y and 6 froma single adjustment of the apparatus, the 
two variables being alternately adjusted until no output is obtained. 
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Fig. 4. Copper films of increasing thickness measured at 45° angle of incidence for various 
wavelengths. The evaporation rate was 1 A per second. 
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Fig. 5. Automatic recording of oxidation of an evaporated film of copper. 


Photoelectric 
Cel 


Bridge to 


Balance Tuned 
Antiphase Amplifier 
; So Signals 
Rotating Auxiliary 
Analyser Analyser 


Fig. 6. 


The principle of the method is illustrated in the schematic diagram fig. 6. 
The general arrangement is similar to that employed in Method A, except that the 
light falling on the photocells is modulated by a rotating analyser placed immedi- 
ately before the double image prism. ‘The modulations of the two beams falling 
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on the photocells contain components of frequencies 2f and 4/, where / is the fre- 
quency of rotation of the rotating analyser. "The 4f components are rejected by a 
filter, and the two photocell outputs thus effectively consist of signals of frequency 
2f. ‘The complete analysis shows that the two signals have, in general, an arbitrary 
phase difference which depends on the mutual orientation of the axes of the vibra- 
tion ellipse, and the axes of the double image prism. When, and only when, the 
axes of the ellipse coincide with the axes of the double image prism are the two signals 
inanti-phase. When this adjustment has been made, the ratio of the two signals 1s 
equal to the intensity ratio of the axes of the ellipse. Thus, by correctly orienting 
the double image prism to coincide with the axes of the ellipse, and balancing the 
signals by means of the bridge or auxiliary analyser, the orientation and ratio of the 
axes of the ellipse can be determined. 


4.1. Sensitivity and Accuracy — 


For comparison with the discussion of Method A in §3.1 above, it will again be 
assumed that an auxiliary analyser is used. A detailed analysis of Method B 
shows that its sensitivity is identical with that of Method A. ‘Thus eqns. (5) and 
(6) apply also to this method. The final accuracy will, however, be affected by the 
additional sources of errors inherent in the use of a rotating analyser which were 
discussed in §2.1. When the method is applied to the direct determination of % 
and A we have a situation analogous to that discussed in §2.1. These special 
effects result in errors in the measured value of %, whereas such errors in A are 
eliminated in a complete set of readings. Similarly, when the method is used to 
determine y and f the resultant errors in the measurement of y are not eliminated 


4.2. Results 


Some results are illustrated in table 3, from measurements on an aluminium 
surface. ‘The apparatus was arranged to measure y and A directly, % being 
obtained from the polarizer adjustment, and A from the setting of the helical 
potentiometers in the bridge circuit. 


Table 3 
80° 00’ 76° 28’ Zoe O00" 70° 00’ 652 007 60° 00’ OO" 50° 00’ 
41° 49’ 41° 20’ = 41° 28’ 41° 44’ 42° 14’ 42° 437 43° 06’ 43° 33’ 
70° 08’ 90" 00°F * 96°"107 114° 42’ 128° 36’ 139° 44’ 148° 16’ NESS G28 


* Values of 6 and y measured using rotating analyser as described in § 2. 


a 


§5. CONCLUSIONS 

The work described in this paper has demonstrated the feasibility of new 
methéds for the analysis of elliptically polarized light. One of the main advantages 
of these methods is that they are applicable to a wide spectral range, provided 
polarizing devices and radiation detectors are available. Moreover, one very 
important advantage which is gained from the elimination of optical compensators 
is that the apparatus does not require recalibration at each wavelength used. 

In our work, the accuracy has been limited mainly by instrumental considera- 
tions (e.g. the divided circles available in our apparatus), and it is not yet known 
whether the extremely high theoretical accuracy suggested by the sensitivity 
calculations can be achieved in practice. On accuracy considerations alone, the 
discussion given above suggests that the use of a rotating analyser is undesirable, 
and the type of system used in Method A seems preferable on practically all counts. 
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It may be advantageous to list here some of the experimental details which require 
careful attention if the highest accuracy is to be attained. 

(i) The linearity of the photocells and their associated circuits are important 
when the ratio is obtained electrically. On the other hand, if an auxiliary analyser 
is used it would appear that only careful matching of the two circuits is required. 

(ii) The circuits, asa whole, must be stable over periods of at least a few minutes, 
during which time y and f may be determined. 

(iii) The light beam must be of uniform intensity and suffer no translational 
movement as the mounting is rotated; alternatively, the photocathodes must be 
of uniform sensitivity. 

The methods described are certainly as rapid as normal visual methods of 
comparable accuracy, and outside the visible spectral range they are much more 
rapid than previous methods. Figure 5 demonstrates the way in which rapid 
changes may be followed by the development of automatic recording devices. 
Such devices may, for example, be used to record changes in A provided that no 
dichroism is present. When the axes of a birefringent material are at 45° to the 
incident light and if a quarter-wave plate is used to reduce the elliptical vibration to 
a linear one, then a small change in A will simply rotate this plane of polarization. 
This change in orientation may be readily recorded by means of the out-of- 
balance signal, which can be shown to be directly proportional to the change in A 
over quite a reasonable range 5A. The method is used to its best advantage 
under these conditions of maximum illumination. The application to accurate 
rapid polarimetry where f is zero and y alone is to be measured is obvious, and 
would enable optical rotations to be determined to greater accuracy than by visual 
means. 
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ABSTRACT. Noble metals show absorption in the blue or ultra-violet; to this is due the 


) colour of copper and gold. As was pointed out by Mott, the main absorption seems due 


to internal photoelectric effect, with excitation of a d electron. This absorption could be 
analysed—as for the alkalis—as a transition between two bands, described by Bloch orbitals, 
here d ands. But we must take into account the correlation between the hole created in the 
nd shell and the Fermi electrons. We can assume that the hole is screened by a (n+-1)s 
electron bound to it, so that absorption should begin for the frequency of formation of the 
excited nd®{n+1)s configuration of the monovalent ion. Agreement with experiment 
obtained using this model is satisfactory both in the pure metals and their alloys. 'The same 
model may be used for x-rays, and optical absorption energy is related to the energy separ- 
ation between the s and d parts of K and Ly; yy emission bands. An order of magnitude 
for the van der Waals interaction-energy between d shells could perhaps be deduced from 
the structure on the short wavelength side of the absorption peak. 


§1. OPTICAL ABSORPTION IN ALKALI AND NOBLE METALS 

HE small optical absorption of the alkali metals is successfully explained 

in terms of the conduction electrons. ‘These electrons may absorb energy 

in two ways (Mott and Jones 1936). First they oscillate in the alternating 
field of the incident photon, and these oscillations are damped by the resistivity 
of the metal. Also, if the energy of the photon is large enough, some conduction 
electrons may make a transition to a state of an upper band in the crystal. ‘The 
two mechanisms seem to explain satisfactorily the optical absorption in alkali 
metals, as was recently shown by Butcher (1951). 

The noble metals on the other hand show a strong absorption in the optical 
or near ultra-violet region (fig. 1, taken from Seitz 1940, after Minor 1903 and 
Meier 1910). ‘This absorption is more than ten times as large as for the alkalis 
and we must look for another mechanism to explain it. 
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Fig. 1. Optical absorption n« in noble metals Fig. 2. Electronic transitions for the 
(after Seitz). edge of the absorption peak in copper. 


Mott and Jones (1936) first pointed out that the three noble metals have 
d shells with not too high ionization energies, and that the strong optical 
absorption was probably related to some excitation of the d shells. 
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Thus in a simple description (fig. 2), we take an electron from a d shell and 
put it at the top of the conduction band. We know from Fuchs’ study (1935) 
that the d shells overlap very little in copper. The energy required to extract — 
a d electron and send it outside the metal (potential 0 in fig. 2) would be much the | 
same as if the Cut ion was im vacuo, but for the fact that the ion is in the (nearly 
constant) field due to the lattice (other Cut ions and Fermi gas). ‘Thus the 
energy required is the difference between the second ionization potential Ey | 
andacoulomb term E,. We do not, however, send the electron outside the metal, | 
but to the top of the conduction band. So we gain the negative energy — Ey 
corresponding to the top of the conduction band. Ey is obtained by the 
approximate relation (Mott and Jones 1936): Ey=E,+#,—0-4 Ey where 
E,, E, and Ey are the sublimation energy, the first ionization potential and the | 
Fermi energy of the metal. 

The total energy required is finally: 

hv=f ES Ey: ON eS eee (1) 
E, is the difference of the coulomb interactions with the rest of the lattice of a 


Cut and a Cu?+ ion. For Cut, contributions by the surrounding atomic | 


polyhedra cancel out very nearly; there remains only the interaction of Cut 
with the part of the Fermi gas enclosed in its own atomic polyhedron. ‘The same 
is true for Cu?+ if we assume that the local density of the Fermi gas remains 
unaltered during the absorption. Computing the fields of Cut and Cu?* ions 
by Fermi’s approximation of the Thomas—Fermi method (Gombas 1949), we 
find then the values of FE, listed in table 1 for copper; those for silver and gold 
are also given. ‘The computed values of the absorption edge Av are somewhat 
larger than the experimental ones. 


Table 1. Optical Absorption Edges Ay (in ev) computed from eqn. (1) 


Metal Ey Ey Ee hvcomp hvexp 
Cu 20°34 8-38 8-04 3-92 2A. 
Ag 21-4 8-33 7:36 5°71 4-0 
Au 19-95 10-87 6:66 2-42 ~2°3 


This could be explained by taking into account the width of the d band. 

As, however, energy £, is large, the ionization of the absorbing ion must 
perturb strongly the Fermi gas, and our picture requires some further refinement. 
This will improve the computed values of the absorption edge and give a possible 
interpretation of its fine structure. 


§2. SCREENING OF THE POSITIVE HOLE 
In fig. 2, we assigned the hole to the d shell of a definite atom. This is not 
quite true: the hole may jump from atom to atom. But it travels very slowly; 
its effective mass is very large or, we may say, the d band is very narrow. The 


hole behaves therefore like a slow and heavy atomic nucleus, in the field of which | 


the conduction electrons, with their effective mass nearly equal to unity, must 


at every instant be in approximate equilibrium.* They must in particular | 


screen the coulomb field of the hole; for it is well known that no such field may _ 


exist in a metal in equilibrium. And the screening most likely is by a bound | 


4s electron (fig. 2). 
* We apply therefore the Franck—Condon principle. 
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We picture therefore the edge of the absorption peak in copper as follows : 
a 3d electron makes a transition to the top of the conduction band, leaving a hole 
in the d shell; the hole is screened by a 4s bound electron provided by the 
conduction band.* 

The 4s electron does not overlap the neighbouring ions much; so its energy 
is not very different from that of the 4s electron of a free Cut (3d° 4s) ion. The 
energy of the conduction electrons on the other hand is probably not very much 
perturbed by the replacement of a normal 3d!° Cut ion by an excited 3d® 4s one. 
The minimum energy required for the absorption must therefore be near to the 
energy of excitation 3d!°+3d%4s of a free Cut ion. 


§3. COMPARISON WITH EXPERIMENT 

Table 2 gives the comparison for copper, silver and gold. Experimental 
values are taken from Seitz (1940), after Minor (1903) and Meier (1910). 
E,, E,, E; and E’ are the nd!°+nd®(n+1)s excitation energies for a free 
monovalent ion, as measured in spark spectra (Cu and Ag, Bacher and Goudsmit 
1932; gold, Platt and Sawyer 1941). The four values correspond to the various 
triplet and singlet configurations of the final state, ?D3, *D,, 3D, and 1D, 
respectively. As we shall explain in a further paper, the excitation to the triplet 
state, which leads to a configuration with a non-zero total spin, is probably weaker 
than that to the singlet state. We should therefore expect an absorption 
increasing rather slowly up to a maximum corresponding to E’. 


Table 2 
Experimental Values (ev) Theoretical Values (ev) 
Edge A Peak B B—A 1B EB, 12 KE E—-E£, 
Cu DA. 2:7 0-5; A Dye DN) Bow (ESS 
Ag 4-0 4-7 0:7 4-83 5:02 5-40 5:70 0:87 
Au ~2°3 3-2; ~1-0 1°85 2:18 3:42 3:66 1-81 


The agreement with experiment of these roughly estimated values is 
surprisingly good. ‘The energy terms neglected here are mainly an attraction 
of the bound state by the neighbouring d shells, and its exchange interaction 
with the conduction electrons, both of which tend to reduce the energy of 
absorption (cf. Friedel 1952). In contrast with what happens in alkali metals, 
both terms should be small, for the 3d*4s ion has a spatial distribution not very 
different from the 3d!° one. This is in agreement with a very rough estimation 
we have made of those terms, following the method indicated in our previous paper. 

Table 2 shows that the difference in energy between the edge and the peak 
of absorption, given theoretically by E’ — Ey, increases regularly with the atomic 
number, as one should expect from a triplet-singlet separation. Gold on the 
other hand absorbs in the optical region, and not in the far ultra-violet as one 
would expect by comparison with copper and silver; this is due to a property 
of the 5d! shell of Aut, probably in connection with the existence of 4f excited 
states of low energy. 

* The 4s bound state may possibly be actually replaced by a heaping up of charge in the 
conduction band, as described by Friedel (1952). ‘This point is not of great importance here; for 


we showed that the heaping up of charge has then nearly the same spatial distribution and average 
energy as the bound state would have if it existed. 
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In conclusion, the main term in the energy of the absorption edge seems to 
correspond to the nd!°> nd°(n + 1)s excitation of the free monovalent ion. This 
is independent of the position of the conduction band. Thus one should expect 
the absorption edge to be very little displaced when alloying elements of valency 
other than unity are added to the noble metal. This seems in agreement with 
experiment (aluminium in copper, McPherson 1940, zinc in copper, Lowery, 
Wilkinson and Smare 1937, nickel in copper, Bor, Hobson and Wood 1939). 
The small displacement of the absorption edge may be due to changes in coulomb 
and exchange interactions of the bound state with the conduction electrons or the 
neighbouring atomic cores. 


§4. DOUBLET STRUCTURES IN X-RAY EMISSION BANDS 

The same sort of description is also applicable to any x-ray transition. At the 
end of a K emission process in copper for instance, the hole in the K shell has 
been filled by an electron coming either from the conduction band (fig. 3 (a)) 
or from the 3d shell (fig. 3(6)); in that case, the hole in the 3d shell must be 
screened as before by a 4s bound state produced by the conduction band. 
Configurations (a) and (6), drawn in fig. 3 for the transitions of highest energy 
(emission edges), have obviously different energies; (a) and (b) correspond 
respectively to the beginning and the end of the optical transition described in 
fig. 2, and their energies differ by that necessary for the optical transition. 
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Fig. 3. End of CuK emission. 
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Fig. 4. Cul? emission (after Bearden and Friedman). Fig. 5. Structure in Cu optical absorption. 


We must therefore expect an emission band with two edges, as is indeed 
observed (fig. 4, from Bearden and Friedman 1940), the difference of energy 
between the two edges being equal to the energy of the edge of the optical 
absorption peak. The passage from the (1s3d!°)Cut ion to the (1s23d9) one is _ 
made possible by the presence of the conduction electrons, in the way described — 
above for optical transitions. Some direct quadrupole transition may however 
occur in the emission of type (6), and explain its rather symmetrical shape. 
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This makes difficult an accurate measurement of the separation between the 
edges. Bearden and Friedman’s data give about 3 ev, which must be compared 
with the 2-1 ev obtained in optical absorption (table 1). 

The Cu Ly, yr emission, where the description is analogous, gives a similar 
doublet, with a separation of about 2:5 to 3 ev (Cauchois 1950, 1952, cf. also Saur 
1936 and Farineau 1939). Here a direct transition of the excited (2p°3d!°)Cu+ 
ion to its stable state 2p*3d° is possible, and gives a very strong (b) peak. 


§5. SCREENING BY p ELECTRONS AND VAN DER WAALS 
INTERACTION BETWEEN d SHELLS 


We have assumed so far a screening of the md hole by a(n + 1)s bound state. But 
screenings by higher excited states are a priori possible, (n+1)p being the most 
favourable energetically. 

We must therefore expect a structure on the short wavelength side of the 
main absorption peak. Meier’s measurements on copper (1910) show a 
secondary peak with a maximum probably near to 2400A or 5:-I5ev. The 
energy for the 3d!°+3d%4p excitation of Cut im vacuo is definitely larger, from 
8:2 to 9-lev (Bacher and Goudsmit 1932). However, the 4p bound state, 
although it has a nodal plane, probably overlaps more than the 4s state does the 
neighbouring copper ions. The resulting attraction by the neighbouring ions 
would explain the difference from 8-2 to 5-lev..- 

If we neglect the absorption of higher frequencies the van der Waals interaction 
energy due to the d shells is given by (Mott and Jones 1936): 

RE S205 38/e*h*\ aif? Zhe 

Seat OEE (=) | ey sf hvhv'(hv+ hv’) | (2) 
Here hy and hy’ are the energies of the 3d'+>3d%4p and3d%4s_ transitions. 
respectively, f and f’ their oscillator strengths and Q the atomic volume. 
3d!°+ 3d%4p is the first direct excitation allowed for the d shell (fig. 5). So the 
first term in (1) gives the interaction between d shells, the second the 
polarization by a d shell of the conduction band on top of a neighbouring d shell. 
The term in f’?/(hv’)? does not appear in (2), for it would correspond to the 
polarization of the conduction electrons by themselves, for which a formula of 
type (2) is certainly not valid and should be replaced by a direct study of 
correlation. 

For copper, hyv=5-l,ev and hv’ =2-7ev; hence FE, ~ —1-7,f(f+5/'). The 
exact evaluation of f and f’ would be very difficult. But we should obtain a very 
rough order of magnitude by considering the absorption peaks as broadened 
lines; f is then given in terms of the maximum absorption u« (Heitler 1936) by: 

3m hy 

t= Gite? N 

where hy is the energy of the absorption peak, y its width and N the number 

of atoms per cm®. From fig. 1, we deduce the reasonable data: nx =nk’ ~2-0 and 
y2-1; y'=1-35ev. Thus f=0-43; f’=0-27, and E, = —1-3,ev. 

This is only a very rough order of magnitude. But its large value may explain 
the discrepancy between the cohesive energy of copper and Fuchs’ computations 
(1935). Adding to Fuchs’ result (Ep;=—1-43ev/atom) a correction due to 
Hsu (1951) for coulomb interaction with the conduction electrons of the slightly 
interpenetrating d shells (Eq = —0-5 ev/atom) we have Ey + Ey = — 1-9, ev/atom, 


ynk 


which differs from the experimental value (Kubaschewski and Evans 
1951) E,=—3-54+0-0,ev/atom by an amount of the order of EY: 
E,—(Ey+Ey)=—1-6ev/atom. A further investigation of the absorption of 
noble metals in the ultra-violet would therefore be of interest, to check our 
interpretation of the structure and eventually to obtain a more precise estimate 
of the van der Waals interactions. 
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The Optical Constants of Thin Metallic Films Deposited by 
Evaporation 


BYAP olen CLEGG 
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ABSTRACT, 'The optical properties of evaporated ARS of silver, gold, tin and indium 
have been determined by a photoelectric method. The influence of factors concerned with 
the evaporation procedure has been investigated in the case of the silver films. ‘The 
experimental results are accounted for by the agglomerated state of the films, and the 


application of Maxwell Garnett’s theory is shown to lead to qualitative agreement with the 
experimental values. 
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§1. INTRODUCTION 
HE measurements detailed below were made with the photoelectric apparatus 
. that has already been described as Method A by Archard, Clegg and 
Taylor (1952). The metallic films were produced by downward 
evaporation on to a clean, wedged, glass substrate at a pressure of about 
2x10->mm of mercury. The required optical properties of the film were 


then determined without breaking the vacuum; these were tan %=7,/r,, the ratio 


of the Fresnel reflection coefficients of light polarized in and perpendicular 
to the plane of incidence; and A=6,—6,, the differential phase change of these 
components on reflection from the film. 

A schematic diagram of the apparatus used is shown in fig. 1 and is largely 
self-explanatory. Measurements were made at angles of incidence of 45° and 65°. 
Without breaking the vacuum this angle could be varied by rotating the 
polarizer and analyser arms, in opposite senses, about the axes DD and BB. 
A 240-watt box-type mercury lamp was used as a source. 
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The mass per square centimetre of the deposited film was directly measured 
by weighing with the vacuum microbalance, described by Clegg and Crook (1952), 
placed adjacent to the substrate. For convenience this measure has been 
transformed into a thickness, the density of the film being assumed equal to 
that of the bulk material; this thickness is known as the ‘ bulk thickness’ of the 
film. For thicker films this value was checked, after their removal from the 
vacuum chamber, by a multiple-beam interferometric method. 


§2. EXPERIMENTAL RESULTS 

The first metal to be extensively studied was silver, chosen because it is inert, 
and is therefore not likely to undergo any chemical change in the atmosphere of 
the vacuum system. Various factors connected with the production of the 
films were altered so as to determine those that were most important. ‘The 
only factors to have any appreciable effect on the resultant optical properties 
were the evaporation rate, the temperature of the substrate and the cleanliness of 
the substrate. 

It was thought that the different topographies of plate glass, fire-polished 
glass and an optical flat might have influenced the growth of the films, but there 
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was no difference between the optical properties of films deposited on various 
types of glass substrate. ‘There was also no appreciable difference between films | 
deposited from tantalum, molybdenum, or plantinized tungsten filaments: it — 
was therefore decided to use the latter for the subsequent depositions. 

In fig. 2 are shown some results obtained for various evaporation conditions — 
at a medium rate of about 2A per second. The scatter of the experimental points — 
for any given run is small, and the scatter between runs gives an indication of © 
the experimental accuracy of these and the other experimental results quoted 
in this paper. 
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An important fact, which enabled measurements to be taken rapidly, was that 
the properties of films of a given thickness were independent of the number of 
stages in which they had been evaporated. For example, the properties of 
films 150A thick built up in 2, 6 or 8 steps were the same. Although this applied 
only to films deposited at, or slower than, an evaporation rate of 2A per second 
it enabled the optical properties of the films built up zz vacuo to be determined 
as functions of film thickness. 

The original measurements on silver films deposited at a rate of about 
2A per second were made at an angle of incidence 6,=65° for a wavelength of 
54614. The results obtained are compared in figs. 3(a) and 3(4) with those 
obtained by Avery (1949); the theoretical curves in these figures were calculated 
from the standard formulae for thin films as given by Drude (1894) and Lord 
Rayleigh (1902), the optical constants of the film (homogeneous) being assumed 
to be those of bulk silver. ‘The disagreement between the author’s results for A 
and those expected theoretically and obtained, down to about 100A, by Avery 
led to further experiments in which the evaporation rate was varied. These 
results are also shown in fig. 3(a), and it may be seen that an increase in the. 
evaporation rate extends the agreement between theory and experiment to 
smaller film thicknesses, though for the thinnest films there is still a major 
discrepancy. 

In fig. 4 are shown for various wavelengths the detailed results for silver at 
6, =65°, the evaporation rate being 2 A per second; it will be noted that there is a 
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complete change in the character of the results, so far as the values of A are 
concerned, between A=5461 A and 43584. A similar change is also apparent in 
the results shown in fig. 5 for an angle of incidence 6, =45° at the same evaporation 
rate; the rise in the value of A above 180° is again quite unexpected; theory, as 
was the case for 6, =65°, predicts the very reverse. From a knowledge of the 
optical constants of almost any bulk metal one would expect for 6, =45° that the 
value of A would fall below 180° as the film thickness increased from zero. For 
convenience one may consider films which do not conform to this expectation, 
ie. those which have a ‘hump’ in their A values, to possess markedly 
‘non-metallic’ optical constants. For both angles of incidence it is only for 
wavelengths of 4358 A and below that theory and experiment even approximately 
agree. 
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The size of the ‘hump’ in, say, \=5461 4 is very sensitive to changes in the 
substrate temperature. The higher the temperature of the substrate during 
evaporation the wider the ‘hump’; the ‘hump’ is also enlarged when the 
deposit is evaporated on to a deliberately contaminated substrate. 

In fig. 6 results for A at 6, =45° are shown for films of gold, tin and indium 
evaporated at rates of about 2, 1-5 and 1-5A per second respectively. It will be 
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noted that the A values for the gold films all fall below 180°, in very approximate 
agreement with theoretical expectations. Both tin and indium, however, show 
large ‘humps’ which extend, especially for the latter, to far greater thicknesses 
than is the case for silver. 
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§3. INTERPRETATION OF THE RESULTS 

The direct determination of the optical constants of a thin film is tedious, 
for, it must be remembered, one cannot derive m and k, where the complex index 
of the film is N=n—zk, explicitly from measured values of #, A and ¢ the film 
thickness, though the converse problem is easily soluble. The method adopted 
by the author was to calculate from the standard equations families of curves, 
with various parameters of n and k, for various values of t. ‘The experimental 
curve was then fitted over a limited thickness range, by extrapolation if necessary, 
to one of the calculated curves, and thus an approximate value for the optical 
constants of the film was found. 

If the results for silver are considered, the most striking fact is the sudden 
change between A=5461 A and 4358A. ‘The two experimental curves for these 
wavelengths are shown, together with some calculated curves, in fig. 7. It will 
be seen that constants 4—1-0z and 2-5 —2-5z, which may be contrasted with the 
approximate values 0-1—3-32 and 0-1—2-372 for bulk silver, give curves lying 
adjacent to the experimental ones over the first few Angstréms, but that a given 
experimental thickness corresponds to a greater thickness on the theoretical 
curve. ‘The implication of this is that the density of the evaporated film is less 
than that of the bulk material. If the factor g denotes the volume of the bulk 
metal per unit volume of the evaporated film, then in this case we see that 
q =0-6, since a thickness of 10A as determined by the microbalance corresponds 
to a theoretical thickness of between 16 and 174. 

The great difference between the values found experimentally for thin films 
and those for the bulk material suggests that the physical state of the metal may 
be different in films and in bulk. Very thin films are agglomerated into discrete 
‘islands’ of metal; this has been shown by the electron microscope studies of 
numerous workers including Picard and Duffendack (1943), Hass (1942) and 
Sennett-and Scott (1950). Maxwell Garnett (1904, 1906) proposed a theory 
originally intended to explain the colours of metal glasses, which he extended 
to deal with the problem of thin metallic films, assumed to be either homogeneously 
mixed with some dielectric material (air or vacuum) or to be composed of spherica 
particles of dimensions smaller than the wavelength of light. 
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Maxwell Garnett’s theory is not, perhaps, the most refined, but owing to 
the indeterminate physical nature of an evaporated film it is considered by the 
author that a simple theory is in many ways the best. Results similar to those of 
| Garnett may be obtained by an extension of the well-known Lorenz—Lorentz 

equation. The additive properties of molecular refractivity defined by this 
equation may also be applied to mixtures if a proportionality factor W is included 
for each constituent 7. For such a mixture we may write 


SW Ly 1) 
et28) nel. ee ( 
where « is the dielectric constant and S is the density. 


For a mixture of two components, both of unit permeability, eqn. (1) 
reduces to 


eg Oe ae Ny? —1 
nad ne@toee aphP Ot ae fio we alee ae (2) 


where 7, is the resulting refractive index and g, and q are the fractional volumes 
of the components a and b. 
Rearranging eqn. (2) and substituting for g,=1—4q,, we see that 


ie a Ny” pat ne ~~, he (3) 
ne+2 qa 242 Ble “at y eT OTM Ran Gawe 


The more complete theory of Maxwell Garnett leads to the following relation : 
ny? = My? Me? oe My 
ne + ony Ya Ny? + 2n," 
which applies to a film composed of spherical particles. 
If one considers a metallic film of index n,=N=n-—t1k deposited in vacuo, 
then ,=1 and eqns. (3) and (4) both reduce to 
2 2 
a =45aG5° Re ide (5) 


which is the form used for the calculations presented in this paper. ‘The theory 
thus enables the resultant optical constants of such an ideally agglomerated 
film to be calculated for various values of g, provided that the optical constants 
of the bulk material are known. Depending on these initial constants, a suitable 
choice of g may well give resultant constants which are markedly different. 

In the case of silver, for example, it can be seen from fig. 8 that for a value 
of g=0-6 the resultant constants of the film change from ‘metallic’ ones into 
‘non-metallic’ ones. Further calculation of the variation of these constants 
with wavelength shows that for g=0-65 the resultant constants of such a film 
would change from about 4:7—0-82 at A=5461A to 2:-5—3-5¢ at A=4358A; 
these values are similar to those found experimentally for such a value of g (fig. 7). 

Further confirmation of the fact that an initial constant of the order 
4— 1-07 is correct for \=5461A is seen from fig. 9, in which the experimental 
values for 6, =45° and 0, =65° are shown together with some selected calculated 
points. It will be noted that as the film thickness increases so the constants 
required to fit simultaneously on both curves tend to the values of the bulk 
constants, i.e. 2 decreases whilst k increases. 

Normally, the anomalous ‘non-metallic’ value of the optical constants for 
silver does not apply to any but the thinnest films, as can be seen from fig. 5. 
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This is presumably because the aggregates rapidly tend to coalesce and the 
constants thus tend towards those of the bulk material. If, for some reason, such 
as the raised temperature of the substrate, the aggregation persists to greater 
thicknesses, then the anomalous ‘ non-metallic’ constants will similarly persist— 
this can be seen to be the case for the silver films deposited on a heated substrate, 
for the ‘hump’ in the A curve is much pronounced (fig. 10). Similarly, it may 
be argued that a slow evaporation rate or deposition on to a contaminated substrate 
favours thé aggregated state and consequently makes the films more non-metallic” 


than normal. 
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By a similar application of Maxwell Garnett’s theory it may be shown that 
the resultant optical constants of an agglomerated gold film for q greater than 0-6, 
although altering largely, do not suffer such a sudden and complete change with 
varying wavelength as do those of silver. For a value of qg equal to 0-6 table 1 


Table 1. ‘The Optical Constants of thin (104A) Gold Films (q¢=0-6) 


A (A) 5790 5461 4358 
Approx. obs. values 2°5—3-01 1-8—2:-52 1-6—1-87 
‘Theoretical values 2°3—3°-91 1:-4—2:47 1:6—1:02 
Bulk gold 0-23 —2-38:7 0-41 —2:-052 1:29—1:677 


shows the experimental results together with those calculated theoretically, and, 
for comparison, the values of bulk gold as found by the author. As the ‘hump’ 
in the A curve for tin is so pronounced, agreement between experimental and 
Maxwell Garnett theoretical values was found up to thicknesses of 254. 
Table 2 shows a comparison of results. 
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Table 2. The Optical Constants of thin (25 A) Tin Films (¢ =0-6) 


A (A) 5461 3560 
Approx. obs, values 3-0—0:57 2:3—1-07 
Theoretical values 2°8—0:32 2:2—1:52 


It will be seen that there is in general good agreement between the 
-experimental results and those predicted on the basis of an agglomerated film; 
as the real film is usually a complex system of aggregates and not very like the 
ideal film considered by Maxwell Garnett, it is somewhat surprising that the 
agreement is as good as it, in fact, is. The character of the results for tin and 
indium, being ‘non-metallic’ to fairly great thicknesses, suggests that these 
metals are aggregated up to thicknesses considerably greater than is the case for 
silver ; indeed, electrical measurements on films of both metals deposited at room 
temperature showed that they were practically non-conducting, at least for 
thicknesses up to 300A. For thin films of tin some recent electron micrographs 
by Scott, McLaughlan and Sennett (1950) show the enhanced nature of the 
agglomerates very well. 


§4.. CONCLUSION 

It has been shown that the application of Maxwell Garnett’s theory of thin 
films is capable of explaining broadly the observed optical properties of evaporated 
metallic films. This implies that by far the most important factor which alters 
the optical constants of thin films from those of the bulk material is their 
aggregated structure. ‘The permanence, as regards film thickness, of such a 
structure may be inferred from the observed optical properties of the film; this 
is made specially obvious when the bulk constants of the metal are such as to give 
“non-metallic’ constants for the thin aggregated films. ‘The evaporation rate 
is found to be of importance. 
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ABSTRACT. When the diameter of the focal spot of electrons on an x-ray target is reduced 
to the order of a micron, the conditions of heat dissipation become so favourable that a greatly 
increased specific loading is possible. The use of an electron lens to obtain such a small 
spot introduces a limitation on the rate of supply of energy, as its high spherical aberration 
requires insertion of a very small aperture, the size of which must be reduced pari passu with 
the size of spot desired. For magnetic lenses of usual design the rate of energy supply is 
then proportional to the radius of the spot. Comparison with the rate of heat dissipation 
shows that, for spots of the order of a micron, the limitation on x-ray intensity is not the 
thermal properties of the target, as in tubes with large focus, but the emission obtainable 
from the cathode. Specific loadings of the order of 10’ watts cm~? can be tolerated in a 
spot of 1 micron; the corresponding emission current density would be of the order of. | 


100 amp cm~?. 


§1. INTRODUCTION 


N the normal type of x-ray tube, with a source of the order of 1 mm? in area, 
| the main limitation on output intensity is the rate of heat dissipation in the 

target. As the size of the spot on the target is reduced the conditions of heat 
flow become more favourable (Miiller 1927, 1929, 1931, Goldsztaub 1947, 
Oosterkamp 1948 a, De Barr and MacArthur 1950). Several types of micro- 
focus tube have been made for purposes where high intensity rather. than large 
total output is important (Ehrenberg and Spear 1951, Witty and Wood 1950). 
Special conditions obtain when a magnetic lens is used to obtain a very small 
spot and when the target is in the form of a Lenard window of thickness comparable 
with the range of the incident electrons, so that heat flow is almost entirely radial. 


§2. HEAT DISSIPATION 


When the radius of the electron spot 7, is large compared with the target thick- 
ness t, heat flow can be regarded as purely axial from the front surface at temperature 
T, to the rear surface at T,. If the thermal conductivity of the target material is 
k, the rate of dissipation of energy is 

E,=7FW,,=27,2(T,— T,)k/t, cles sn 
where W,, is the intensity or specific loading. When 7, is of the same order as 
t radial heat flow becomes important. Oosterkamp (1948 b) obtains the combined 
energy flow EF, by introducing into eqn. (1) a correction factor s, dependent on the 
ratio of y, tot. For a constant value of s, as obtained by imposing the condition 
ry =t, we have E,cr, and W,1/r,. ‘The same simple relation holds when r, <t. 

When the target is so thin that electrons penetrate through it, as in a micro- 
focus tube recently constructed (Cosslett and Nixon 1952), energy is supplied to 
the metal throughout an approximately cylindrical volume. Assuming that this 
supply is uniform along the path ¢, and that the front and rear faces of the target are 
at the same temperature 7}, the flow of heat will be purely radial and 


EB, =k Thea T\antin(Rin he ee ce 
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where R is the radius at which the target has the temperature 7. In practice 
convection at the rear face of the target will lower its temperature below T, ; 
radiation in all practical conditions is negligible. 

The figure shows the variation with spot diameter of E, and Ep for a copper 
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Rate of supply and dissipation of energy for a target spot of diameter d, and a 10 kv electron beam. 
Es=rate of supply for given value of 8B. ER=rate of dissipation by radial heat flow for given 
target radius R. He=rate of dissipation by combined radial and axial flow. 


target of thickness 1p. It is assumed that 7, =100°c and 7, =800°c; the same 
results will obtain for a tungsten target at T,=2000°c. ‘The two curves for Ep 
assume values of R=0-05 and 0-5 mm respectively. 


§3. HEAT SUPPLY 

As in optics, the energy flux in an electron lens system is constant across any 
plane normal to the axis, that is, 8, the energy per unit area per unit solid angle, is 
invariant with axial distance for a given beam voltage and in the absence of 
absorption and scattering. If p is the current density at any plane at which the 
potential is ¢ and the angular width of the beama, thenB=p/7«?. Ifr,isthe radius 
of the image, the power incident on it is 

Eas Pxhe = 1 Ta top Os =CONSt.,, ne wane (3) 
SINCE 79% is invariant, according to the electron optical form of the Helmholtz— 
Lagrange equation. Thus, if an electron lens were perfect, the rate of energy 
supply would remain the same whatever the size of the image spot formed, whereas 
the rate of dissipation in a metal target would fall with spot size according to eqn. (1) 
or eqn. (2), depending on the conditions. 

In fact magnetic electron lenses suffer from severe spherical aberration, for 
which no correction is known. ‘The aperture of the imaging beam has to be 
limited to make the radius of the disc of confusion 6, less than that of the desired 
spot. Ifalens has aspherical aberration coefficient C,, then 


pss EC ge od ih re (4) 
In practice it is sufficient to prescribe r, =6,, as the disc of minimum confusion is 
smaller than that in the gaussian image plane. Substituting in eqn. (3), 
te eae eee Dy (eS wate y| ie |y Wy ligerty sor sbre ote (5) 


784 V. E. Cosslett 


C,, which is of linear dimensions and of the order of the focal length f, is constant 
only at fixed power in a magnetic lens ; as the excitation is changed, C, varies 
rapidly. In the case of some simple field distributions it is found that C, oc f*, but 
for actual lenses it follows from recent measurements (Liebmann and Grad 1951) 
that C,ocf?? over a wide range of high power. Equation (5) then yields E,« 72; 
in the figure E, is given for two values of specific emission f at a beam voltage of 
10 ky. 
§4. CONDITIONS IN THE TARGET 


The thermal conditions in the target will depend on whether the lens is operated 
with fixed aperture or with x, always adjusted to satisfy eqn. (4) for the chosen spot 
size. In the latter case we have E,ocr, over a range of lens power practically 
down to the minimum focal length. If both radial and axial flow contribute to 
the dissipation of heat eqn. (1) with Oosterkamp’s correction gives E,<7, also. 
Hence target temperature 7, is independent of 7,. The value of 6 which would 
raise the target to a temperature 7, is found by setting E, equal to £, and inserting ~ 
the experimental dimensions and beam voltage. ‘The conditions which give £, as 
in the figure demand for a 1 spot an emission 6 = 6-7 x 10° amp cm *steradian 1, 
which is much greater than the specific emission of a tungsten filament in normal 
operating conditions ; to raise the target to its maximum permissible temperature 
would require an incident intensity of 7-5 x 10?ampcm™. For this to be supplied 
through the aperture fixed by eqn. (4) the source would have to emit 30 amp cm”. 
By overrunning the filament the emission could be raised towards the value needed 
to make optimum use of the thermal properties of the target, but at the cost of 
reducing the filament life to the order of a few hours. Increase in accelerating 
voltage ¢,, if permitted by other considerations, might be a preferable alternative 
(eqn. (3)). 

Conditions become worse after the minimum focal length of the lens has been 
reached, at given ¢,. Further reduction in spot size can only be obtained by 
reducing the size of the electron source or increasing its distance from the lens. The 
rate of energy supply is then given by eqn. (5) with C, constant, so that E,«7,5* 
and diverges rapidly from the limiting rate of dissipation F, if 8 is constant. 

In practice a lens will often be operated with an aperture of fixed size, such that 
the spherical aberration is less than the minimum spot size required. Then 
in eqn. (3) a, is fixed and E,ocr,”. The rate of energy supply falls faster, with 
diminishing spot size, than does the rate of dissipation which, from eqn. (1) with 
Oosterkamp’s correction, follows E,ocr,. For a given maximum 7; a condition 
of energy balance will be reached at a spot size depending on the limiting value of f. 
The tube will be most efficient in x-ray production if 8 can be made great enough 
to raise the target spot to the maximum permissible temperature at the smallest 
value of r, required. 

As shown above, this is possible below 27, =1,: only at the cost of shortened 
filament life. If it can be realized, then E£, will be much larger than £, 
for appreciably larger target spots, and unless the filament emission is reduced 
when the spot size is increased the target may reach the melting point. 


§5. CONCLUSIONS 


A target spot smaller than about 1p can only be raised to the maximum safe 
operating temperature by so increasing the intensity of the electron source that 
its lifetime is appreciably reduced. As the spot size is reduced a limit will be 
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reached in which the lifetime of the filament becomes equal to the required exposure. 

The emissive power of the filament is thus the limiting factor on x-ray intensity and 
not the thermal properties of the target, as in a normal tube. The criterion in 
choosing a target metal then becomes efficiency of x-ray production, instead of 
thermal conductivity and melting point. Metals of high atomic number, such as 
tungsten and gold, which also have good thermal constants, are then to be preferred 
to copper. 
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The Optical Properties of Liquid Selenium 
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Communicated by R. W. Ditchburn; MS. received 8th April 1952, and in final form 
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ABSTRACT. ‘The optical properties of liquid selenium in the near infra-red have been 
measured for a temperature range up to about 400° c. ‘The absorption edge moves out to 
longer wavelengths as the temperature increases by an amount corresponding to a change 
of energy gap of 1:4x10-%ev/c°. Absorption is very small at low frequencies. For 
long wavelengths the refractive index decreases as the temperature increases by an amount 
corresponding to the thermal expansion. 


selenium falls rapidly with increasing wavelength in the red end of the visible 

spectrum and is very low in the infra-red. ‘The photo-conductivity threshold 
does not coincide with the absorption edge (Weimer 1950, Weimer and Cope 1951, 
Gilleo 1951). Both the absorption edge and the photo-conductivity threshold 
move towards shorter wavelengths with decreasing temperature. As far as is 
known, the absorption of liquid selenium has not previously been measured. 
Since the structures of liquid and amorphous selenium are believed to be similar 
in nature, it is to be expected that liquid selenium will show optical properties very 
similar to those of amorphous selenium. 

The absorption of liquid selenium as a function of wavelength was measured for 
a temperature range up to about 400°c. ‘The selenium was contained in fused 
silica absorption cells which fitted closely into an electrically heated brass block. 
The thicknesses of cell used were 5mm, 1cmand2cm. ‘The cells were placed in 
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|: has been shown by many workers that the optical absorption of amorphous 
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a parallel beam of light at the entrance slit of an infra-red monochromator. The 
light source was a tungsten filament lamp, the light from which was interrupted at a 
frequency of 800c/s. The receptor cell was a lead sulphide photo-conductive 
cell which was coupled to a high-gain tuned amplifier. Figure 1 shows trans- 
mission against wavelength for a cell of 2cm thickness. ‘The movement of the 


absorption edge towards longer wavelengths as the temperature increases is very _ 


marked: in addition the edge (plotted in terms of wavelength) becomes less sharp. 
The absorption coefficients were obtained by the usual method of plotting 
log (transmission) against cell thickness. Figure 2 shows absorption coefficients 
plotted against frequency. The curve retains much the same shape over the whole 
of the temperature range. Using the present technique it was not possible to 
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measure absorption coefficients greater than about 20cm. ‘The fact that absorp- 
tion is very small at low frequencies indicates that the free carrier concentration is 
small; this is confirmed by the fact that liquid selenium possesses a high thermo- 
power (Borelius and Gullberg 1944). 

The absorption edge must be related to the energy difference between two states 
of the material, but since the edge does not coincide with the photo-conductivity 
threshold this energy difference does not correspond in a perfectly direct way with 
a ‘forbidden gap’ width. Defining an energy E49 =A where vo is the frequency 
corresponding to an absorption coefficient of 20cm™!, it is found that 
Ey9 =2°16—1-4 x 10-8T ev (T in °K). Ina similar way E,.,=1-86 —1-4 x 1037. 
Two determinations of the photo-conductivity threshold by Gilleo (1951) satisfy 
the equation Ep =2-62 — 1:4 x 10-37. 

The refractive index of liquid selenium was measured by using a cell with 
faces inclined at an angle of about 20 degrees. The cell mounted in the heating 
block was placed on the table of a spectrometer of which the collimating slit 
coincided with the exit slit of the monochromator. ‘The eyepiece of the spectro- 


meter was replaced by the lead sulphide cell. ‘The refractive index as a function of | 


temperature and wavelength is shown in fig. 3. The values for the amorphous 


; 4 
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‘Wd solid at 18°c agree to within 1% with previous measurements (Dowd 1951, 


Gebbie and Saker 1951). When the index is known the amount of radiation lost 
by reflection from the two silica-selenium and the two silica-air surfaces of the 
absorption cell can be calculated. The result of this calculation agrees to within a 
few per cent with a measurement of the amount of reflected radiation obtained by 
extrapolating to zero thickness a graph of log (transmission) against thickness. 
For long wavelengths where the dispersion is small the variation of [7] with temper- 
ature (where [7] =(n?—1)/(n?+2) and m is the refractive index) is such that 
(1/[n]) d[n]/dT =1-7 x 10-*per°c. The writer has measured the expansion of 
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Fig. 3. The refractive index of liquid selenium. 


solid and liquid amorphous selenium by interference and dilatometer methods 
respectively, and for both types of material the volume expansion coefficient is 
1-5x10-*per°c. The refractive index therefore appears to depend on the 
density p in the manner expected from theoretical considerations, so that 
(n* — 1)/(n? + 2)p is constant with temperature. 
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The Contamination in Evaporated Films by the Material 
of the Source 


By O. S. HEAVENS 
Department of Physics, University of Reading 


Communicated by R. W. Ditchburn; MS. received 23rd March 1952 


ABSTRACT. Two methods have been used to detect contamination in evaporated films } 
of silver and germanium by the boats or filaments used as sources. It is shown that a | 
limit of a few parts in ten thousand of tungsten or molybdenum is attainable by microchemical 
methods. ‘The much lower limit of one or two parts in 10’ is attained by the radioactive | 
tracer method. The extent of contamination of films of silver and germanium by tungsten | 
and tantalum boats has been studied using the tracer technique and the minimum contam-_ 


ination attainable found to be of the order of a few parts per million. The possibilities of | 
transfer of the boat material by direct evaporation or by reaction with the residual gas in the | 


‘system have been minimized; the results suggest that contamination arises by the solution 
of the boat materials by the molten silver or germanium. 


§1. INTRODUCTION 
q HE increasing importance of films produced by the thermal evaporation 
process makes it desirable to know the extent of contamination of the films 
by the material of the boat or filament under the normal evaporation 
conditions. [wo methods of determining this contamination have been 
studied: (i) microchemical analysis of the evaporate, (ii) use of radioactive tracer 
source elements. As would be expected, the latter method is considerably more 
sensitive than the former and is the only method capable of detecting contamin- 
ations smaller than 1 part in 10®. The microchemical method can nevertheless 
provide a qualitative indication of the presence of impurities to the extent of a 
few parts intenthousand. For cases in which this limit of detection is sufficiently 
low, the chemical method is to be preferred on account of its speed and simplicity. 
The methods have been used to study the purity of films of silver and 
germanium when prepared by evaporation from boats or filaments of tungsten, 
molybdenum and tantalum. ‘Two types of vacuum system have been used for 
these experiments. For tests at pressures down to 5<10-*mmHg, a metal 
system employing an oil diffusion pump was used; for pressures down to 
10-? mm Hg, a glass system with a mercury pump and solid CO, + acetone was 
used, 


§2. MICROCHEMICAL TESTS 
Using the methods given by Feigl (1947) films of silver prepared by evaporation 
(a) from a molybdenum boat and (4) from a tungsten filament were examined for. 
contamination. ‘The boats or filaments were first thoroughly cleaned and out- 
gassed in vacuum at 2000°c and were then loaded with ‘Specpure’ silver. 
Films of area 2 cm? and thickness 3 000-5 000 A were prepared at three evaporation 
temperatures, measured by optical pyrometer. ‘The films were removed from the 


glass slide by the minimum necessary amount of aqua regia, on which the spot 
tests were made. 


Determination of Purity of Evaporated Films 789 


The thickness of the films was measured by placing a glass slide, partly covered 


| with a razor blade, beside the target surface, thus enabling a portion of film with a 


well-defined edge to be collected. This slide was subsequently covered with an 
opaque silver layer and the step height determined using multiple-beam Fizeau 
fringes (‘Tolansky 1948). This method has been shown to be reliable for films. 


_ of this order of thickness (Heavens 1951). The electron micrographs of Sennet 


and Scott (1950) show that the density of films of this thickness is unlikely to. 


| differ-appreciably from that of the bulk material. 


For the detection of molybdenum use is made of the brick-red coloration 


») which results from the addition of stannous chloride to a mixture of a molybdate 


solution and potassium thiocyanate. A test solution containing a known amount. 
of molybdenum showed that the identification limit (0-ly) and concentration 
limit (1: 500000) quoted by Feigl could be reproduced under the conditions of 
these tests. The corresponding test for tungsten is much less sensitive, the 
detection limit being 5y. 

An alternative test makes use of the catalysing effect of tungstates on the 
decolorizing reaction of titanium trichloride with malachite green. ‘This action 
is markedly accelerated by the presence of small amounts of tungstates in acid 
solution. From test solutions containing known concentrations of tungsten, 
weighed drops containing 0-ly to 0-5y of tungsten were added to microdrops of 
1°, titanium trichloride and 0-005°%, malachite green on a spot plate. In the: 
absence of any added tungstate the green coloration was found to disappear in 
from 3 to 10 minutes: drops containing 0-2—-0-3y of tungsten caused decoloration 
in from 1 to2 minutes. Results for drops containing 0-1y were erratic, indicating” 
0-2-0-3y as the reliable detectable limit. 

No tungsten or molybdenum were detected in silver films for which the: 
evaporation temperature was less than 1700°c. For the temperature range 
1750-1 800° c, just detectable amounts of the boat material were observed with. 
somefilms. Taking 0-2yas the detection limit in each case the observed contamin- 
ation at 1 750—1800°c corresponds to 2-3 parts in 104. This figure is confirmed. 
by the results of the radioactive tracer tests described below. 


§3. RADIOACTIVE TRACER METHOD 

Many properties of evaporated films are likely to be only negligibly affected by 
smaller amounts of contamination than are detectable by the microchemical 
methods discussed above. It is known, however, that the electrical properties of 
semiconducting materials are very sensitive to minute amounts of impurity, 
concentrations far below the limit detectable by chemical methods being sufficient 
to cause significant changes. For certain types of impurity, the use of the radio- 
active tracer technique enables the detection limit to be lowered by two or three 
orders of magnitude below the chemical limit. The availability of the method 
depends on the existence of an active isotope which has a suitable half-life and 
which can be produced in sufficient concentration. Measurements have been 
made using tungsten and tantalum, for which elements very suitable active 
isotopes may be produced by (n, y) reactions. ‘The tantalum isotope 1®?'T'a has a 
half-life of 117 days: since the whole of natural tantalum consists of 1*"'T'a, a large 
activity (80 mc/g) is obtainable. The detection limit attained with this material 
is less than 10-°g. With the tungsten isotope 1®’W, of which the half-life is 
24-1 hours, an activity of 360 mc/g is obtained, yielding a detection limit of less. 
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than 10-1°g. For molybdenum, the tracer method offers but little advantage over » 
the chemical method. The isotope **Mo possesses a suitable half-life (67 hours) 
but the low activity at present attainable enables a detection limit of only 10-*g to | 
be reached. 

Small boats, made from 0-2 mm sheet, of tungsten and tantalum were irradiated 
in the Harwell pile until sufficient activity had resulted. ‘The boats were then 
used as evaporation sources and the evaporated materials tested for active elements 
in the usual way. 


§4. EXPERIMENTAL DETAILS 
(i) The Metal System | 
The metal system is a standard vacuum evaporation plant with Pyrex bell-jar. | 
By means of suitable long-handled tools, the active boat is installed and clamped. — 
A platinum—platinum-rhodium thermocouple is mounted on a light spring so 
that when the boat is installed, the couple is held against the bottom of the boat. | 
(The activity of the boat makes the usual spot-welding process of attachment too | 
hazardous.) An optical pyrometer is used for temperatures above 1600°c. 
A shutter is moved over the boat to collect the first part of the evaporate which may | 
consist of surface contamination, for example from oxidation by the atmosphere. 
All vacuum taps, shutters etc. are handled by remote control. 


(ii) The Glass System 
The glass system consists of a large U-tube mounted horizontally (tig. 1). 
‘The tungsten boat is mounted on the end of two tungsten rods, 3 mm diameter and 
60 cm long, which are prevented from sagging by glass bridges. The targets 


Iron Slug 


— —— 


to Pump 
Fig. 1 


are nickel discs of the type used in the standard Geiger counter assembly. — In the 
glass system these are held on a glass spoon which is kept out of sight of the boat 
during preliminary outgassing. A magnet enables the target to be moved into 
position to receive the evaporated material. 


Determination of the sensitivity. 

A minute (~100y) weighed piece of boat material is attached to aluminium > 
foil with collodion and irradiated in the pile simultaneously with the boat. An 
identical aluminium foil plus collodion is also irradiated so that the (small) activity 
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due to these supporting materials may be readily allowed for. (‘This activity was 
found to decay to negligible proportions in a few hours.) The use of a small 
weighed specimen, for which the disintegration rate is within reach of the normal 
counting apparatus enables the specific activity of the boats to be found easily. 

The sensitivity obtained with the tantalum was such that N counts per minute 
above background corresponded to 2:3,x10-°x Ng. The slow decay of the 
Ta isotope allows the fall in sensitivity over a period of a few days to be 
neglected. With tungsten the initial sensitivity was such that N counts per minute 
corresponded to 9-5 x 10-2 x Ng. 


$5. RESULTS OF TESTS 


(i) Tests with Boats Alone 


A preliminary series of tests was made in which the boats, after initial cleaning 
and outgassing for several minutes at over 2000°c, were heated empty to various 
temperatures from 1 000° cto 1 900° cand the targets tested for radioactive material. 
The amounts of tungsten and tantalum transferred to the target were found to be 
considerably larger than would be expected from direct evaporation at the temper- 
atures used. The mass of metal transferred was found to depend directly on the 
pressure in the system, results obtained for tungsten at 1200°c being shown in 
fig. 2. Substantially similar results were obtained with the tantalum boat, the 
rate of loss of tantalum ranging from 5 x 10-§gcm™~*sec“ at 5 x 10-*mm Hg, to 
45 x 10-$gcm- sect at 4x10-°>mmHg. Over the range of pressures used, the 
rate of loss of boat material is directly proportional to pressure, suggesting that the 
mechanism of transfer involves reaction with the residual gas in the system. 
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Fig. 2 


The two most likely mechanisms would seem to be (i) oxidation of the boat by 
residual oxygen, (ii) oxidation by reaction with residual water vapour. ‘The oxides 
of tungsten and tantalum thus produced evaporate immediately at temperatures 
above 1000°c. Langmuir (1915) has studied the first-mentioned process in the 
case of tungsten and finds that only a fraction of the incident oxygen atoms reacts 
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with the tungsten. At a pressure of 10-?>mm and temperature of 1200°c, | 


Langmuir’s figures lead to a rate of loss of tungsten of 5 x 10-'gcm*sec™ for 
this process, compared with the observed value of 1-2 x 10-8 gcm™?sec”” if it be 
assumed that the residual gas is air. 

It is thus clear that reaction with residual oxygen will not account for the 
observed transfer of boat material. Moreover, observations on the nature of 
residual gases in vacuum systems of this type show that the predominant constitu- 
ent is usually water vapour (Blears 1951). ‘Tungsten is known to react readily 
with water vapour at the temperatures of these experiments (Langmuir 1913) 
so that some transfer of material by this process is to be expected. Ifit be assumed 
that the residual gas is entirely water vapour, then the number of water molecules 
striking the boat per second is sufficient to account for the observed transfer. 
From the point of view of studying the contamination of germanium and silver 
films, it is seen that, provided the pressure in the system is kept below 10-mm, 
no direct transfer of tungsten occurs within the limits of detection at an evaporation 
temperature. of 1200°c. Any tungsten detected in films produced at these 
pressures must therefore be carried up with the evaporated material, and further 
reduction of pressure in the system is unlikely to produce films of lower tungsten 
content. 

(it) Tests with Silver and Germanium 


Weighed quantities of these materials were evaporated completely from the’ 


active boats and the mass of the films produced determined either directly by 
weighing or interferometrically. ‘The silver or germanium charge was first 
melted so that the whole of the surface of the boat was covered. During the last 
stages of the evaporation, part of the boat surface became exposed. The time for 
which any of the heated part of the boat was exposed to the target did not exceed 
about 5°% of the total time of evaporation, due to the effective wetting of the boat 
materials by the molten silver and germanium. Blank experiments using empty 
boats were performed under the same conditions of temperature and pressure to 
determine the mass of boat material transferred to the target by the process 
described in (i) above. At all temperatures, the mass of boat material in the silver 
and germanium films was found to be considerably higher than the mass transferred 
by reaction with residual gas in the system, under the same conditions. 

The results of the tests are given in tables 1 and 2. 

The figures in parentheses in the last columns show the amount of contamin- 
ation which would be produced by the mass of boat material transferred by reaction 
with the residual gas of the system, as indicated by the blank experiments. (These 
figures represent an upper limit since most of the boat surface is covered by the 
evaporating material during the test.) The much larger observed values of 
contamination suggest that solution of the boat materials by the molten silver and 
germanium occurs even at the lowest practicable evaporation temperatures. 
At higher temperatures, the degree of contamination increases considerably ; 


silver films containing over 10% of tantalum have been produced by evaporating — 


at temperatures of about 2500°c. 

In some experiments at the lower temperatures (1000°c-1200°c), the 
evaporation was stopped before any of the boat surface became exposed. The 
figures obtained were found not to differ significantly from those in which the boat 
contents were evaporated completely and have been included in the mean values 
for each temperature which are given in tables 1 and 2. 
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Table 1. Tests in Metal System (pressure 8-15 x 10-* mm) 
Boat Material Evap. Metal Evap.Temp.(°c) Contamination of Film x 108 


Ww Ag 1000 4-9+0:3 (<0°5) 
1100 22 (1) 
1410 120 (20) 
1900 330 (90) 
Ge 1070 23403 (<0°5) 
1250 5+ (1-5) 
1470 64 (27) 
1800 430 (105) 
Ta Ag 1100 20 (3) 
1200 29 (3) 
1410 45 (10) 
1670 113 (29) 
Ge 1120 11 (3) 
1380 yi (4) 
1530 17 (12) 
1800 76 (37) 


Table 2. Tests in Glass System (pressure below 5 x 10-7 mm Hg) 
Tungsten Boat 


Evap. Metal Evap. Temp. (°c) | Contamination of Film x 105 

Ag 1080 (Peep (0) 
1200 34 (0) 
1510 68 (6) 

1740 106 (19) 
Ge 1040 2°14+0°2 (0) 
1200 4+3 (0) 

1620 39 (10) 

: 1800 330 (42) 


It appears therefore that the smallest degree of contamination attainable in 
films of silver and germanium evaporated from tungsten or tantalum boats is of 
the order of a few parts per million. If it is desired to study properties of films 
which are likely to be affected by contamination of this order, then an alternative 
method of evaporation suggested by Ashworth (private communication) in which 
the heating is effected by radio frequency currents may be used. 


REFERENCES 


Buears, J., 1951, 7. Sci. Instrum, Suppl. No. 1, 36. 

| Fgic., F., 1947, Qualitative Analysis by Spot Tests (New York : Elsevier). 

Heavens, O. S., 1951, Proc. Phys. Soc. B, 64, 419. 

Lanemuir, I., 1913, Trans. Instn. Elect. Engrs., 32, 1893 ; 1915, J. Amer. Chem. Soc., 37, 
417. 

SENNET, R. S., and Scott, G. D., 1950, 7. Opt. Soc. Amer., 40, 203. 

TOLansky, S., 1948, Multiple Beam Interferometry (Oxford : University Press). 


794 


Self-Absorption in Arc Sources in Thermal Equilibrium 


By: H: -EDELS 
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ABSTRACT. An analytical determination of the self-absorption of a spectral line emitted 
by a source is not normally possible if account is taken of the spatial variations in the source 
of the excited atom concentrations and spectral line profiles. The problem is simpler for 
uniformly excited sources, for example sources in thermal equilibrium, but not for high- 
pressure arc sources normally assumed to be in thermal equilibrium, because of the existence 
of temperature gradients. In the paper an analysis is developed which shows that if such 
arcs are treated as uniformly excited, an error of less than 5°% is introduced into the deter- 
mination of the self-absorption of spectral lines, emitted by transitions between an upper 
state s and a lower absorbing state t, which satisfy the limits hvgq/RT, >10 and vg/r4q <0°5, 
where q is the ground state and J, the axial source temperature. The theory is then applied 
to the experimental data for a high-pressure mercury arc discharge and to new measurements 
of the line widths at 5461 and 4358A. It is shown that the calculated value of the 
self-absorption at 5461 A agrees with the experimental value only when the line has a 
resonance profile. The existence of this form of profile is also indicated by considerations 
of the broadening processes in the discharge. 


§1. INTRODUCTION 
HEN radiation of a given frequency is emitted by any small volume of 
a gaseous light source it has normally to pass through a portion of the 
source before being received by the observer. During this process a 


fraction of the emitted radiation energy is absorbed by atoms and molecules: 


similar to those causing the emission. ‘This phenomenon is termed self-absorption, 
and a knowledge of its magnitude is often required in spectral investigations of 
gaseous electrical discharges. 

In order to determine the self-absorption account must be taken of the spatial 
variations of excited atom concentrations and spectral line profiles. The 
complexity of the problem is such that analytical determinations are not 
normally possible, and recourse must be made to graphical methods. Such 
n-ethods are, however, laborious, so that simplified source conditions are normally 
assumed in order to obtain analytic solutions. ‘This is evident in the comprehensive 
work of Cowan and Dieke (1948), who obtain the general equation and consider 
its application to the characteristics of different source models and the extent to 
which these models are comparable with a practical source. In the exhaustive 
treatment of Bartels (1949, 1950), although a rigorous solution of the problem 
is not achieved, an approximate general solution is obtained which is applied 
to the special case of self-reversal, when the reversal peaks have reached their 
maximum value. 

From the work of these and previous investigators it is clear that the 
mathematical difficulties arise from the condition that the concentration and 
line profile variations of the emitting and absorbing states are different, and that 
the problem is considerably simplified for uniformly excited sources, that is, 
sources throughout which the ratio of the concentrations of emitting and 
absorbing atoms is constant. ‘This condition holds, for example, for sources in 
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thermal equilibrium, in which the energy of all types of particle are distributed 
according to the Maxwell—Boltzmann law. However, in arc sources which are 
accepted as being in a state of thermal equilibrium, for example in the carbon 
are in air, as shown by Ornstein and Brinkman (1934) and Kruithof and Smit 
(1944), and the mercury arc by Elenbaas (1951) (see also Edels 1950), appreciable 
temperature gradients exist, as shown, for example, by Kruithof (1943), Fischer 
and K®6nig (1938), Elenbaas (1951), so that the assumption of thermal equilibrium 
is invalid when applied to the whole discharge region. The following analysis 
of self-absorption shows that, despite the existence of temperature gradients, it 
is often possible to assume that the source is uniformly excited without causing 
undue errors in the final results. This makes possible self-absorption calculations 
of the more complex cases of these sources, including line profile variation. 
Such calculations are of considerable importance because the presence of 
self-absorption in a source in thermal equilibrium produces emergent radiation 
which is not descriptive of a thermal source. 


$2. DERIVATION OF A SIMPLIFIED SELF-ABSORPTION EQUATION 

We will first obtain the general equation for self-absorption and derive from 
it the simplified equation which holds for the special case of sources in thermal 
equilibrium. 

Consider a light source in a steady state so that the emission and absorption 
processes have no effect on the concentrations of the atoms in the excited states. 
If we then take a small volume dv of the source containing dN, atoms/cm? in 
an upper state s capable of emitting radiation by transition to a lower state t with 
a frequency between v and v+dy, and dN, atoms/cm? in the lower state capable 
of absorbing this radiation, and if this volume is subject to an intensity of 
radiation /,, then the energy received from the solid angle dw is 

A(E,dv) = {dN,Ay—(dN,B,—4N,By)1,}hv dwdv/4m ...... (1) 
where Ax, (ByJ,) and (B,,/,) are Einstein’s probabilities of spontaneous and 
induced emission and absorption respectively. Now in most investigations the 
solid angle of the observed radiation from the source is small, and the radiation 
may be taken to exist as a plane wave directed along the axis of observation x. 
‘Thus the change of intensity over dx is 

hv d 2 d 

d(I, dv) = {aN Ang - Beet (1 aM ae) 1} fc ae ca (2) 
where the relationships between the Einstein coefficients have been substituted 
and g, and g are the statistical weights of the excited states. In most practical 
cases Kirchhoff’s law is obeyed, so that the spectral absorption and emission 


frequency distributions are the same at any point. ‘Thus, assuming a symmetric 
source about x =0 such that (N,),=(N4g)of(x) and (N;)2=(Ni)of (x), we have 


(dN,)o=(Noofl)flv, x) dv and (dNi)e=(NdoF(a)f04%) dry cee. (3). 


where (NV,)) and (N;)p are the total concentrations of excited atoms in the states s 
and t at x =0, and a normalization condition exists given by | f(v, x) dv =1 over the 
range — co to +00. Substituting (3) into eqn. (2) we obtain the general 
self-absorption equation 


d(I, dv) ={(I,)oflx)f(v, x) dv —KyF(x)f(v, x)(I,dv)}dx ecu. (4) 
where (Le)o a (N3)oA ¢/tv/47, Ko = {(M)oAstr*e,/8794} {1 a (N;/+) a&t/Es}- 
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(I,)p is the true intensity emitted at « =0 and K, is the integral of the absorption 
coefficient over the whole spectral line at x =0, as given, for example, by Mitchell 
and Zemansky (1934). The term 1— (N,N, abulde is in most practical cases 
nearly unity and may be ignored, although with the approximation which is to 
be made it may be reintroduced later. 

If the variations along the x axis are assumed to be even functions and we 
take a boundary for the line of frequency v of +X, then the solution of eqn. (4) 
gives the received intensity in the ewe interval dv as 


Ts dv=(ldudefesp (—KuideD f°” L fexp (Kop) rex [K(2GH)°—¥)} i‘ 
(5) 


where (5)”=|*F(x)f(v, x) dx. With the general sbinte in this form it will be 
apparent that the difficulties in applying it to practical cases are due almost 
entirely to the presence of the factor {f(«)/F(x)} and that the expression is 
appreciably simplified for a uniformly excited source since f(x)/F(x) is a constant. 
If we now consider a source which has a temperature gradient of such a magnitude 
that the source still approximates to a state of thermal equilibrium, then it has 
been shown (Edels 1951) that {F(«)}?=f(x) if hy, q/RT)>10 and vg¢/r4q<0-5 
where 8 =1,,/v4g, Where q represents the ground state and 7) is the temperature 
at the source centre. The limit of vy/%, precludes the application of this — 
relationship to resonance lines, a restriction which appears also in Bartels’ work. 
With {F(«)}°=f(x) in eqn. (5) the solutions are still laborious. However, for . 
many spectral lines 6 does not differ greatly from unity, thus the Balmer series 
gives Bnay~1-3, and since it is the integrated values of F(x) and f(x) which affect 
absorption calculations, it would appear feasible that no serious error will result 
by putting 6=1. That such an approximation is possible is shown by the 
following calculations. 

To simplify the calculations assume that the spectral line profile does not 
vary spatially and consider only that half of the source furthest from the observer, 
with «=O at the furthest point. Thus with {F(«)}*=f(«) we have 


T dual gdh Weep | LK fie) | : F(x) as | : {F(«)¥8 exp | Ki fon F(x) as | dx. 


If we now take a function F(x) =(«/X)" and assume a linear spectral profile of 
half-width Av and central frequency vp) given by f(v) =(1/Av){1 — |(v—v9)/Av] }, 
then from (6) the total received intensity J, as a fraction of the true emitted 
intensity Z, which would be received if no absorption were taking place, is 

given by ; 


Ge Suet eine b8 . 
en f bs 243 3(G+2)~ PLAST) | HHLANSE GLA) +H] 


where 6=XK,/Av(n+1) and 6=n(B—1)/(n+1). When B=1, 4=0, and the 


fraction of the intensity received is given by 


if b ihe 68 
0-7-[1-247-aten-- Rod OG (8) 


The maximum error in assuming 8=1 for a given value of } is obtained 
when ¢ is a maximum; for hydrogen this occurs when 8 = 1-3, so that with n =4, 
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¢=0-1. In this case when b=1, Q=0-736 and P=0-746, giving an error of 
slightly greater than 1°% when the absorption taking place is of the order of 26%. 
When 6=2 the absorption is 43° and the error is only 2:3°%%. Now for the 
maximum value of v,4/v,q4 assumed previously, the corresponding maximum for 8 
is 1-5. With this value and b= 1 and 2 we obtain P=0-753 and 0-593 respectively, 
giving errors less than 5°%, which are within the range of accuracy normally 
obtained in intensity measurements. ‘Thus for a source in thermal equilibrium 
satisfying the conditions hy,,/RT) > 10 and vg./r4q <0-5, self-absorption calculations 
may be made assuming that F(x) =/(«). With this assumption eqn. (5) reduces to 


Tdy= om —expi—2K (by tev. da ven (9) 


This is the well-known amie for a source of constant excitation having 
a uniform emission (J,)), uniform coefficient Ky and an equivalent absorbing 
depth of 2Av)(b)>*. With this equivalent depth we may utilize relationships 
previously obtained for constant excitation sources. ‘Thus for a source of 
circular cross section of radius R, if we neglect spatial variation of spectral profile 
and consider the absorption along the source diameter we obtain the ratio of the 
received to true emitted intensity J/J, for different profiles as follows: 


Resonance : f)= ~ / E + la (v— oh | ; 


Then I/I,=exp(—ap)[J(tap)—tJ,(tap)] naw (10) 
(Ladenburg and Reiche 1913) where a=(2/z) [i F(z) dz; z=r/R; p=RK,/Av 
and J), J, are Bessel functions of the first kind. 


2 {in2\i2 Zc 2 
Doppler : f= nen (~*) exp| —In2{— (—»)} fi 


ope ap (ap)? 
Then 7=[1- Bie aan | erste (11) 
(Ladenburg 1930) where a=4{(In2)/7}!? [j F(z)dz and p=RK,/Av. 
1 = 
Linear: Jo)= iN E - ee ] , 
bP ag 1 
Then I, =| 1- ap (lexp(—e)} Hae (12) 


where a=2 {j F(z) dz and p=RK,/Avp, by integration of (9). In each case the 
absorption is dependent upon p, the absorption parameter, and in general, for 
any profile varying spatially only in half-width, it can be shown that the absorption 
is always a function of p. ‘Thus, expanding and integrating (9) we obtain the 
total received intensity of a spectral line from vp to vg+ dv as 


I=4da [ {b- Fees See... bay 


The first term of this expression is the true ida intensity J,, so that 


ree 2K, fy? dv ae 2K) J tf, Bento ash ae] 
I, za lm xb dv Pius wie S! PP reyes) ah : 
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Now in general f(v, x) =(c/Av,)f[(v—v)/Av,], and Ay, in terms of Avy the half-— 
width at x=0 is given by Av,=Av)f(x/R). Since F(x) can also be stated as a- 
function. of «/R, substitution in (13) and change of variables gives for the general 
term [xs_,(—RK,/Avp)%~*] where « is a constant. Thus (13) becomes 


RK, RKo\? 
1 [1 2, Re +04 (4) ~... | =01-A0)1 


The absorption parameter, after substitution of the expression for Ky, and 

giving the state concentrations their thermal values, is given by 
P=(RUN 4) AstA*Zs/87g Avo] [1 — exp (— hrg/RTo) |. 

From this it will be seen that an increase in absorption will take place with 
increases in the absorbing state concentration, source dimensions, transition 
probability and wavelength, but decreases with increase in half-width. 

It is of interest to note that the analysis precludes self-reversal, since 
differentiation of eqn. (9) shows that the maximum of the received intensity 
distribution still occurs at the central-frequency. 


§3. SELF-ABSORPTION IN ARC SOURCES 


The simplified self-absorption equation (9) may be applied only to those 
sources in which thermal equilibrium exists and to which the given limits apply, 
so that F(x) may be taken to equal f(x). 

Since most high-pressure arc columns (~ 1 atm) are believed to be in thermal 
equilibrium, it should be possible therefore to obtain further confirmation of the 
validity of the simplified equation from the experimental data for this type of 
source. ‘The data available are, however, extremely meagre: only the carbon 
and mercury arcs have been shown to be in thermal equilibrium, and self- 
absorption data have been obtained only for the latter case. We will use the data 
provided by Elenbaas (1951) for his mercury arc discharge (0-88 atm, 6 amp, 
Tya3 = 5 885° K at x =0), in which it is clear that the lines at 5461 and 4358A are 
self-absorbed... Dr. Elenbaas has kindly measured the half-width values of these 
lines and obtains 0-25A for 5461 A and approximately the same for 4358 A, so 
that only the profile shapes remain to be determined. A good indication of these 
may be obtained from a consideration of the processes affecting the line 
broadening in the discharge, and from the work of other investigators. 

The natural width of the lines may be neglected as well as the Doppler width, 
which is approximately 0-02A at T~5900°K. Since neither Wendt and 
Wetzel (1916) nor Haslam (1935) obtained measurable Stark displacements 
of these lines, we assume that the quadratic Stark effect constants are less than 
10-4. With this value and an electron concentration of 5 x 10!%/cm? obtained 
from Saha’s equation with a reduced ionization potential of 9-7 ev (Schulz 1948), 
we obtain, using Lindholm’s (1945) formula, a negligible Stark width of the 
order of 0-024. It would appear therefore that the lines are mainly pressure 
broadened, and this is confirmed by the work of Rompe and Schulz (1938) and 
Schulz (1938) on mercury arcs at 20 atm. ‘These investigators find that 
asymmetric pressure broadening and resonance profiles are associated with the 
series (n?S,—2°P,), which includes the line at 5461 A. It may be concluded there- 
fore that the lines at 5461 and 4358A from a mercury arc at a pressure of 
approximately 0-9 atm also have a resonance profile. In this case the absorption 
of the lines may be determined from eqn. (10) if the limits of the analysis are 
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satistied. Calculation shows that for 5461 A, Avyq/RT) ~15 and vg./r4q ~ 0-42, 
whilst for 4358A the latter factor is 0-59, exceeding the limit of 0-5. Thus our 
analysis may only be correctly applied to 5461A. In determining K, for this 
line we use the value of A,, as given by Schouten and Smit (1943) reduced by 
one-sixth as suggested by Elenbaas (1951, p. 37), thus A,4,~6 x 107 sec-!, and 
since R=2-05 cm and Av=2-52 x 10! sec, we find p=39. From the radial 
temperature distribution (Elenbaas 1951, p. 73) it is possible to calculate F(z) 
along the discharge diameter for both the emitting and the absorbing states. 
Since the main approximation made in the analysis is that F(z) is the same for 
both states, use of either function in eqn. (10) should produce a negligible change 
in J/J, if the approximation is valid. Calculation gives {§ F(z) dz for the emitting 
and absorbing states as 0-247 and 0-276 respectively, yielding corresponding 
values of J/Z, from eqn. (10) of 0-306 and 0-292. ‘The insensitivity of the 
absorption to F(z) shown here is independent of the spectral profile. ‘Thus 
similar calculations for a Doppler profile using eqn. (11) give the Z/J, values as 
0-113 and 0-103, whilst a linear profile gives 0-10 and 0-09. It is clear from these 
calculations that the absorption is indeed insensitive to the change in F(z) between 
the emitting and absorbing states, thus indicating the validity of the approximation 
made. 

It is evident, however, that the magnitude of the absorption is appreciably 
affected by the form of the spectral profile. It is thus possible to estimate the 
profile at 5461 A by comparing the theoretical values of absorption obtained for 
different profiles with the absorption experimentally determined by Elenbaas 
(1951, p. 145). In his experiment Elenbaas uses two discharges L, and L, and 
a monochromator arranged to measure only the radiation emitted by the discharge 
diameters. By measuring the intensities of radiation from L, and L, separately 
and together, denoted by J, and J, and J, . respectively, a measure of the absorption 
along the diameter of L, of the radiation received from the diameter of L, is given 
by «={1,+1,—J,,}/;. Since sufficient data are available only for the Hg 
discharge (0-88 atm, 6 amp) we will consider only the case when L, and L, are 
similar discharges. 

If the radiation received from L, is (1,dv), given by eqn. (9), then after 
passage through L, it is reduced to the value (J,dv), exp[—2K(x)9"], so that 
the amount absorbed is (J, dv), [1 — exp {—2K,())"}]._ Thus 


r+ 0 +0 
[tex (2K mre [ [7 texp {4K (Wola | 


+0 eS F2ars . 
[ot -exp{-2K(udy | | [1 -exp{—-2K(p).t}] av 


This is the equivalent of the line absorption of Ladenburg and Reiche (1913). 
For a resonance profile we obtain from (14) where a= (2/7) [5 F(z) dz 
he E _ €2t IS ((Ziap) — tJ ,(2tap)} 
yr eS (tap) —1J (ap)} J 
As before, p=39 and with either 0-247 or 0-276 for J)! F(z) dz we obtain from 
this equation «=0-57. This compares favourably with the experimentally 
determined value of 0-54. 
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Closer agreement than this is not to be expected because of uncertainties in 
the transition probability value, the actual profile form and half-width and 


because of the neglect of spatial variation of profile and unknown accuracy of — 


the intensity measurements. These uncertainties make reasonable a value 
of p ~ 18, which for a resonance profile gives full agreement between experiment 
and theory. The calculated values for « for other profiles are widely different 
from the experimental value of 0-54; thus a Doppler form gives 0-90 whilst 
a linear form gives 0:97. The calculations thus indicate that the line 5461 A has 


a profile corresponding to a resonance form, a conclusion in agreement with — 


that obtained from considerations of the broadening processes in the discharge. 
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ABSTRACT. ‘The calculation of ring shims of the kind suggested by Rose has been 
extended to include polecaps in the form of truncated cones. Design data are presented in 
graphical form. Actual and predicted performances are compared for one magnet. 


certain physical phenomena. An example is the study of nuclear 

magnetic resonance absorption, which usually calls for a field of several 
kilogauss, homogeneous within at least 1 part in 10* over a volume of 1 cm’. 
It is customary to meet these requirements by means of a magnet equipped with 
plane parallel soft-iron polecaps provided with ring shims of the kind calculated 
by Rose (1938). 


VV homogeneous magnetic fields are required for the investigation of 


(b) 


Fig. 1. Diametrical sections of typical magnet polecaps. Fig. 2. Conformal mapping of the ¢-plane 
on the ¢-plane. The real ¢ axis 
coincides with the surfaces of the 
polecaps.. 


| 

. Rose’s analysis envisages polecaps having axial symmetry, a diametrical 
| section of which is shown in fig. 1(a). The raised rim round the edges of the 
| polecaps compensates to a first order for the radial fall in field towards the edges 
) of the gap. Rose gives in graphical form the height of the necessary raised rim 
_ in terms of its width. His data apply only to the case shown, namely that for 
which the angle PQR is 90°. 

When a permanent magnet is employed, the efficient use of the permanent 
magnet material frequently requires the diameter of the poles to be greater than 
the diameter of the gap. The polecaps must then be in the form of truncated 
cones (fig. 1(5)). The angle P’Q’R’ is now greater than 90°, so that Rose’s data 
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are not strictly applicable. The present paper gives the results of an extension 
of the calculation to polecaps of general angle. 

It will not be necessary to repeat here Rose’s analysis in so far as it is general. 
It will suffice to outline the method and state his results. He makes two 
assumptions: (a) that the reluctance of the soft iron of which the polecaps are 
made is negligible: the surfaces of the polecaps may then be treated as — 
equipotential surfaces; (b) that the diameter of the polecap is large compared | 
with the width of the gap: the potential problem then approximates to one in | 
two dimensions, and may be solved by conformal mapping using the Schwarz 
transformation method (see for example Kottler 1927 or Green 1948). 

Cartesian coordinates are chosen as shown in fig. 2. The x axis lies in the 
median plane of the gap, and the z axis is the normal to the median plane 
through the edges of the polecaps. Half the width of the gap is taken as the 
unit of length. A complex variable €=«*+7z and an auxiliary variable ¢ are 
introduced, the relationship between which is given by 


dt Y {2 \%/% —-1 |) 
G=-su(1- 4) eed es (1) | 


where the ¢,, are the values of ¢ (all real) at the corners of the polygon formed by 
the air gap, the «, are the internal angles of the polygon, and the product is 
taken over the corners for which ¢t,>0. Three values of ¢ may be disposed 
arbitrarily (say 0, 1, 00) as shown in fig. 2. The values —1 and — © follow from 
the symmetry of the system about the median plane of the gap. Values of 
+t, and +f, are assigned to the remaining corners as shown. ‘The diagram of 
fig. 2 omits the corner P’ of fig. 1(6). The neglect of this corner is justified 
provided that P’Q’ is greater than about half the width of the gap, a condition 
which is usually met in practice. We shall return to this point later. 

Rose shows that the necessary condition for obtaining a first-order correction 
for the inhomogeneity of the field due to the edge of the gap is that the coefficient 
of ¢? shall be zero in the power series expansion of the expression. 


{2 \1-epn/= 
¢.2 
n n 


in which, as before, the product is taken for values of ¢,, greater than 0. This 
requires that 


n 


We will now apply this requirement to the case shown in fig. 2. The angle 8 
is the semi-angle of the truncated cones which form the polecaps, and will be 
expressed as 7. ‘lhe width and thickness of the rim are a and b respectively. 
The parameters of the corners are therefore: 


ig Le OG 

i(1 —b) 1 (3/2 —7) 
a+i(1—b) t, 37/2 
att ty WT are 


Substituting these values of ¢,, in (2) we find 


t= t7 1 + (L—2nlt Allee a eee (3) 
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The value of a is obtained by integrating (1) from 1 to ¢,. This gives 


2/1 pz 1/2 / #2 ~1/2 
= ly eee he )x cat (ee ee ES 
a= = jit (1-22) (5 1) (Fa 1) dt Seahit Ladlad (4) 


The value of 6 is obtained by integrating (1) fromt,toz,. This gives 


aie #2\1/2 / #2 ~1/2 
aes 1/1 _ 42\3—7 Pina 4a 52 
b= = { acre (1. A) (5 1) dials COR (5) 


For a given » and for each value of ¢, between 0 and 1, values of a and 6 are 
obtainable from (3), (4) and (5). Eliminating the auxiliary variable ¢, the 
thickness of the rim is found in terms of its width a. 

The integrals in (4) and (5) have been computed for semi-angles 0 equal 
to 0°, 15°, 30°, 45° and 60° (7 equal to 0, 4, 4, 4, 4). For each @, values of a and b 
were obtained for ¢, equal to 0-4, 0-5, 0-6, 0-8 and 1:0. In fig. 3 the thickness 
of the rim 6 is plotted as a function of its width a for the five values of 6. 
Figures 4 and 5 show respectively curves of a against 6 for constant 5 and of 6 
against @ for constant a, the values being obtained from fig. 3 by interpolation. 


Fig. 3. Dimensions of ring shim. The width of the ring a is plotted against its thickness b for 
various 0, the semi-angle of the truncated cone. The unit of length is half the width of the 
gap. The broken lines are contours of the parameter f,. = 


The integrals were evaluated using Simpson’s rule. ‘The integrand in (4) 
is finite throughout the range of integration. In (5) however the integrand 
becomes infinite at the lower limit ¢,. Simpson’s rule was therefore applied 
from t, down to fy, a value slightly greater than ¢,. For the remaining range of 
integration from fy to ¢, all the factors in the integrand except the last were 
treated as constant, ¢ taking the value ¢,, appropriate to the mid-point of this 
short range. ‘The integral could then be evaluated in this range giving the 
following contribution to b: 


(2t [at ,)(1 — ty)? "(1 — tpn? /t,?)2? Inn [t/t + (to?/to2—-1)¥7]. «..... (6) 
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When a magnet is designed the value of b is often determined by consideration 
of the clearance required in the gap. If the designer wishes to calculate his 
own shim dimensions using (4) and (5) rather than by interpolating from fig. 3, _ 
it is necessary for him to find by trial the value of ¢, which gives the desired value _ 
of b, and then to use this value of ¢, to find the value of a. In order to facilitate _ 
the finding of the value of t, in such cases, lines of constant t, have been drawn — 
in fig. 3 (broken lines). 

We must now return to a consideration of the effect of neglecting the 
corner P’ of fig. 1(6). The surface P’S’ beyond the corner may not be the 
surface of soft iron of low reluctance, but may for instance be permanent magnet 


O16 = r T T T | ae eT T a 


04 


03 0-12 


01 


Fig. 4. The variation of a, the width of the Fig. 5. The variation of b, the thickness of 
ring shim, with the angle 6 for various the ring shim, with the angle @ for various 
thicknesses b of the ring. widths a of the ring. 


material instead. Nevertheless, it should be a satisfactory approximation for 
our purpose to consider that P’S’ has the same magnetic potential as the polecap. 
A value ¢, is assigned to this corner. Reference to (1) shows that an extra factor 
(1 —7?/t,2)"” now appears in the integrand of (4) and (5). An additional term also 
appears in (2), so that (3) must be replaced by 


tat 1 — 291 —t ee (7) 
Our use of eqns. (3), (4) and (5) is equivalent to the assumption that t,-?<1. 


The value of ¢; is found by integrating (1) from 1 to ft, to give the value of the 
length of the conical surface P’Q’ 


2, ts 2 1/2 / 72 —1/2 t2\" 
= igo So | eee wig ate ath yrs oe 
re = | “e1) (5 1) (5 1) (1 =) ditt, ee (8) 


Approximate evaluation of (8) for various 7 shows that for r equal to 1-5 gap 
units t;~7, the actual value depending somewhat on 7 and ¢,. The effect of the 
terms in ¢; * is in fact negligible for r>1-5. For r=1 it is found that tz~4, and 
the terms in ¢;-* become appreciable. For a fixed value of a, the correction to 
be applied to 6 increases with the angle 0; at @=0° the correction is of course 
zero, while at @=60° 6 should be about 5% larger than the value given by fig. 3, 
an increase of the order of 0-003 half-gap unit. The error caused by treating the 
problem as two-dimensional may often be of this order. 

Thus, as stated earlier, provided r exceeds about one half-gap unit the data 
given in figs. 3, 4 and 5 should be satisfactory. If 7 should be less than this, the 
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problem must be worked out specially. The value of t, which gives the correct 
value of 7 for the particular @ must be found by trial from (7) and (8). The 
values of a and 6 may then be found from (4) and (5) after inclusion of the 
additional factor (1 —??/t,”)” in both integrals, using the new value of ¢, from (7). 

A ring shim of the type described here has been incorporated in a large 
permanent magnet designed by the present authors. A brief description of this 
magnet has been given by Hadfield and Mawson (1952), who were responsible 
for its construction. ‘The soft-iron polecaps are in the form of truncated cones 
of semi-angle 38° 40’, having a thickness of 2-5 cm and faces of diameter 20 cm. 
The magnet was designed to have a gap of 5-40cm. The thickness of the ring 
shim was chosen to be 2mm, and for this thickness the appropriate width was 
found to be 5-43mm. After construction the magnet was found to have a gap 
of 5-54cm instead of the specified value of 5-40cm; the dimensions of the shim 
should therefore have been 2-6 % larger. 
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Fig. 6. Comparison of actual and predicted performances for a magnet having 6=38° 40’, and a 
poleface diameter to gap width ratio of 3-6. The relative variation of the magnetic field in 
the median plane is plotted against distance x (in half-gap units) from the edge of the gap. 
Circles: experimental points. Upper curve: predicted variation with shim Lower curve: 
predicted variation without shim. 


The field in the median plane has been measured by search coil and 
fluxmeter and in fig. 6 is shown plotted against x, the distance inwards from 
the edge of the gap. ‘The variation expected theoretically both with and without 
the ring shim has been calculated in the manner described by Rose and is also 
shown in fig. 6. While the agreement with the upper curve is not perfect, it is 
clear that the first order of inhomogeneity has been very largely removed. ‘This 
is particularly gratifying, since it might have been doubted whether the ratio 
of polecap diameter to gap-width, 3-6, was sufficient to satisfy assumption (d) 
adequately ; the relatively small value of this ratio may in fact be the main cause 
of the discrepancy between the experimental points and the upper curve in fig. 6. 


REFERENCES 


GREEN, S. L., 1948, Theory and Use of the Complex Variable (London : Pitman). 
HaprFiELp, D., and Mawson, D. L., 1952, Brit. 7. Appl. Phys., 3, 199. 

Kortter, F., 1927, Handbuch der Physik, 12, 480 (Berlin : Springer). 

Rose, M. E., 1938, Phys. Rev., 53, 715. 


806 


Observations on Growth and Etch Phenomena 
on Haematite (Fe,0;) Crystals 


By A. R. VERMA 
Royal Holloway College, University of London 


Communicated by S. Tolansky ; MS. received 2nd April 1952, and amended 29th May 1952 


ABSTRACT. New observations of growth and etch features on haematite (Fe,O3;) 
crystals are reported. Of the many growth patterns observed, the simplest are the 
elementary ‘ molecular growth spirals’ in accordance with the theory of Burton, Cabrera 
and Frank. The step height between the successive arms of an elementary spiral has been 
measured by the application of multiple beam interference fringes and is shown to be equal 
to the x-ray value of the repeat distance. More complex growth patterns are illustrated and 
discussed. The etch figures reported are also of molecular thicknesses and consist of a very 
large number (~107/cm? of the crystal surface) of similarly oriented small triangular figures, 
and some more complex shapes. 


$1. INTRODUCTION 

CCORDING to the theory of crystal growth developed by Burton, Cabrera_ 

A Frank (1951), when a crystal grows due to the presence of dislocations 

the steps exposed on the crystal face by screw dislocations wind themselves 

into spiral forms giving rise to ‘molecular growth pyramids’. ‘The observation 

of numerous ‘growth spirals’ on the faces of silicon carbide crystals has been 

reported earlier and shown to be in accordance with theory (Verma 1951). Similar 

observations of ‘ growth spirals’ on haematite (Fe,O3) crystals were reported in a 

preliminary note (Verma 1952b), and in the present paper a fuller account of the 
observations of growth and etch phenomena is given. 


§2. TECHNIQUE OF OBSERVATION 
Most of the observed features are of molecular height and are consequently 
faint. ‘The technique found most suitable for their observation is that previously 
employed for silicon carbide crystals (Verma 1951) and consists of depositing a 
highly reflecting layer of silver on the crystal surface by thermal evaporation and 
then examining it in reflection using phase-contrast microscopy. 


§3. STRUCTURE OF HAEMATITE AND ITS MODE OF GROWTH 

The structure of haematite may be referred to rhombohedral axes making an 
angle of 85°42’ with each other or to hexagonal axes with an axial ratio c/a = 1-36. 
The structure assigned to haematite by x-ray measurements (Bragg 1924, Pauling 
and Hendricks 1925) may be pictured as follows. A series of rhombic cells is 
formed by slightly extending cubic cells along a trigonal axis so that the angle 
between the edges meeting in this axis is 85° 42’ instead of 90°. At every cell corner 
is placed a unit consisting of two Fe atoms and three O atoms. The two Fe atoms 
are separated by a distance of about 2-88 A and lie at equal distances from the cell 
corner along the trigonal axis which passes through it. The three oxygen atoms 
surround each cell corner at a distance of about 1:47 A and lie in a plane through the 
corner at right angles to the trigonal axis, i.e. in the (111) plane. The equilateral 
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triangles formed by the three O atoms are oriented at 60° to each other in successive 
layers, i.e. are A and V in orientation. The distance between successive oxygen 
layers is nearly 2-34. Figure 11* shows the structure projected on the (111) 
plane. ‘The structure repeats itself along the trigonal axis after six such (111) 
layers, giving the repeat distance equal to nearly 144. The rhombic unit cell is 
‘outlined by joining the Fe pairs ‘1’ to ‘3’; ‘3’ to{5’; ‘5’ to‘7’; ‘7’ coin- 
-ciding in projection with ‘1’. This structure gives an approximate hexagonal 
‘close packing for oxygen atoms. 
Haematite, which is formed in various ways and is found in rocks of all ages, is 
considered to occur in connection with volcanic activities as a sublimation product, 
often as small thin crystal plates which are parallel to the c axis (Dana 1948). 
At a temperature of about 2 000 °c, at which ferric oxide can be sublimed in crystal- 
line form, it becomes dissociated into a magnetic oxide and consequently it will be 
the magnetic rather than the ferric oxide that will volatilize at these temperatures. 
Haematite is probably produced by the action of water vapour upon ferric chloride 
with the generation of HCl. The HCl vapour, if not removed, can act strongly 
on the haematite crystals formed, and most of the surfaces will easily etch. 


§4. OBSERVATIONS ON GROWTH FEATURES ON HAEMATITE 
Simple or Elementary Type of Spiral 

On the haematite crystals many growth patterns have been observed, some of 
which are of the simple type. Figure 1 shows such a spiral with elementary 
Burgers’ vector (see $5) observed on the basal plane of a thin hexagonal plate of 
haematite. Originating from a right-handed screw dislocation, the spiral is seen 
to terminate on a left-handed one after about two turns. ‘Thereafter successive 
closed loops are generated in accordance with theory. 

This spiral with straight edges indicates that the rate of growth depends 
strongly upon the crystallographic directions, the straight edges marking the 
directions of minimum growth. According to Burton, Cabrera and Frank’s 
theory (1951) the conditions of growth for a polygonal spiral are such that the 
“Frenkel kinks’ or exchange sites are few and far between, whereas the distance 
moved by the adsorbed molecule is small. 

The trigonal symmetry of the spiral in fig. 1 follows from fig. 11, since for 
spirals on the basal planes (111) we are looking on the edges of ‘monomolecular 
layers’ along the triad axis. Further, the spiral (fig. 1) is observed to be triangular 
at the centre, becoming hexagonal after one or two turns. This is represented in 
fig. 12, which is a projection on the (111) plane, remembering that the successive 
arms of the spiral lie on successive ‘ molecular layers’. 

Figure 2 is another example, of which the central growth pyramids are shown 
enlarged in fig. 3. Here two right-handed screw dislocations close to one another 
each give rise to a growth pyramid with very close spacing, indicating that the 
growth of these pyramids took place at quite a high supersaturation. (That these 
features are elevations, i.e. pyramids, has been confirmed by the application of 
multiple-beam fringes of equal chromatic order.) ‘These two step lines join 
together after about one turn (fig. 2) to form a common spiral. Furthermore, the 
step lines are broad and kinked and the corners concave. ‘This is to be attributed 
to the subsequent etching of the crystal giving rise to a general retreat of the step 
lines and rounding off the corners. 

* All the figures in this paper are plates, printed at the end of the issue. 
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Interaction of Growth Fronts originating from Elementary Dislocations 


A striking example of the interaction of the growth fronts from a number of | 
dislocations emerging on the face of a haematite crystal is illustrated in fig. 4. . 
This interaction contrasts with similar observations on silicon carbide crystais. | 
Among the points to be specially noticed are: (i) the meeting of the three edges at : 
an angle of 120° to each other at the centre; (ii) the way in which close pairs of ' 
dislocations give rise to peninsula-like plateaus at two places round the central 
part ; (iii) in a total of about 20 dislocations, only two or three are of the right-handed | 


type, the remainder being left-handed; (iv) by going round the compound growth | 


pattern once we shall descend or move up by as many steps as the effective 
number of dislocations enclosed within the circuit. 


Curvilinear and Bunched Spirals 


On one crystal face a large number of growth hills were observed to be arrayed 


} 


along the edge of the crystal. These growth hills were bunched and their spiral 
form could not be resolved clearly. If these growth hills have originated from 


| 


screw dislocations this would be the simplest case of an assembly of dislocations, — 
and would then form a complex transition surface (Burgers 1937, 1940). On 


the rest of the crystal face similar spiral features (see fig. 5) have been observed. | 


These growth spirals appear to originate from multiple dislocations and might 
have been expected to take on hexagonal shapes with straight edges if left undis- 
turbed. Further, in these spirals on the three alternate edges the steps dissociate 
into several components, giving an appearance analogous to that found in similar 
observations on silicon carbide crystals (Verma 1952a, fig. 3). ‘This observation, 
together with the curvilinear and bunched spirals of high visibility (see fig. 6), 
shows that, like Si-C crystals, it is also possible in haematite for spirals to 
originate from dislocations of multiple strength, which suggests that small portions 
of haematite crystals may exist as polytypes. 

In a few cases rounded spirals which are heavily bunched have also been 
observed. Figure 6 is an example in which the spiral starts at the centre as a 
single step line of high visibility. Over small segments of the spiral the step lines 
have nearly straight edges, especially near the origin ; further from the origin the 
steps dissociate into several components. In the lower half of the figure over a 
dozen elementary faint step lines with the edges straight bunch together, leading 
to curved step lines of high visibility. This suggests that the high visibility 
rounded spiral of fig. 6 may be considered to originate from a dislocation of multiple 
strength. 

In another case a very large number of dislocations (over 100), all of the same 
hand, originate in one region of the crystal. Groups of these step lines bunch 
together, leading again to steps with high visibility. 

The circular type of spiral could be formed when the growth rate does not 
depend upon the crystallographic directions. Bunching of steps seems to 
facilitate the generation of curvilinear spirals. 


§5. MEASUREMENT OF STEP HEIGHT 
Multiple-beam interference fringes (Tolansky 1948) have been utilized for the ' 
measurement of step heights. ‘The individual steps are rather close, and step 
heights are too small to be detected by the shift of the Fizeau fringes at each step. 
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However, as in previous work (Verma 1951), measurements were obtained by 
allowing a Fizeau fringe under high dispersion to pass over the peak of the pyramid, 
and, from the angle through which the fringe appears to bend at the peak, the average 
step height can be determined. Consider a pyramid formed of a series of parallel 
steps with straight edges, having a constant spacing d and step-height h: if a 
Fizeau fringe passes perpendicular to the series of steps on one side of the peak, 
it will appear to bend through an angle 6 such that 


d sin 6 1a 
(1+ 3 cos? @)1/2 X 2 


where X is the distance between the successive Fizeau fringes for wavelength A. 
For small values of @ (which usually correspond to small values of hf) as a first 
approximation this reduces to h=4dsin06(A/2X), which is identical with the 
relation used for circular spirals. 

This method was used for the measurement of the step height of the spiral in 
fig. 1. Since the separation between the successive steps in fig. 1 gradually 
decreases on going away from the centre, the Fizeau fringes passing on these 
steps were observed to be slightly curved. For the evaluation of step height in this 
case the Fizeau fringe picture and the microphotograph at the same magnification 
were superposed, so that the exact position of the fringe was known with certainty. 
For measuring @ two straight portions of the Fizeau fringe, one on either side of the 
centre of the spiral, were chosen and their corresponding spacing of steps d was 
used. 

Using this method, the step height measurements for different independent 
trials were 12:4, 16-3, 11-2, 12-9, 17:2, 13-3, 16-04. 

Another measurement of average step height is possible with the fringes of 
equal chromatic order formed when the slit of the spectrograph was adjusted to 
pass over the centre of the image of the pyramid (fig. 3). This determination 
gave an average value of nearly 11-5 A step height, but is subject to the error that if 
the slit does not pass exactly over the centre of the pyramid the value of step height 
will appear to be low. 

These measurements, within experimental errors, show the step height to be 
equal to the x-ray value of the repeat unit along the trigonal axis. ‘This spiral, 
therefore, is an elementary growth spiral, originating from a dislocation of unit 
Burgers’ vector. 


re 


§6. MOLECULAR ETCH -FIGURES 

Very often on haematite crystal surfaces a variety of etch figures is observable. 
Figure 7 is an enlargement of one part showing a very large number of similarly 
oriented triangular figures covering the whole crystal surface. (For this observa- 
tion the crystal surface was not silvered.) Sometimes, due to the overlapping of 
several smaller triangles, a much bigger triangular figure can be formed. Also 
some hexagonal figures, with a stepped structure and sometimes with a dot at the 
centre, can be observed in fig. 7. This observation of discrete triangular figures at 
~ different points is in accordance with the general observation that solution does not 
act uniformly over the whole face, but begins only at isolated points, proceeding 
more easily in certain directions, giving rise to regular figures. In the present 
case the triangular figure indicates three preferred directions of solution inclined 
to each other at 120°. The triangular shape of the etch figure conforms to the 
trigonal symmetry of the crystal face. 


810 A. R. Verma 


If etch begins favourably at certain isolated points in an otherwise continuous 
lattice surface it would seem that there are small scale lattice irregularities at these 
points, since solution takes place favourably at ridges, steps, etc. The small 
cracks of possibly near molecular dimension assumed by Smekal can explain the 
appearance of etch pits at such spots. According to the ideas of Frank, disloca- 
tions which provide steps for the growth of the crystal should also promote the 
development of etch pits. If, indeed, these etch pits have originated on the ends 
of dislocations, the surface density of these pits would give the density of disloca- 
tions. In fig. 7 this density is of the order of 107/cm?, which may be compared 
with the largest density of dislocations of the order of 105/cm? observed on silicon 
carbide crystals (Verma 1952 a). 

These triangular figures are the elementary etch figures with shallow depths, 
possibly of near molecular depth, a conclusion suggested by the fact that these 
features are not normally observable with usual bright field illumination, and some 
of the triangular features have a low visibility, even under phase-contrast illumin- 
ation. However, because of the smallness of these features, it is not possible to 
apply interferometric techniques for their study. A hollow, which may be a 
different type of etch pit, is shown in fig. 8. This consists of a series of roughly 
circular steps some of which are seen to interlace with one another. ‘The step 
height is approximately 100A. Fringes of equal chromatic order over this 
feature show it to be a pit whose total depth is about 1 000 4. 

Figure 9 is a further illustration ; here an anticlockwise triangular spiral is 
surrounded by the rugged and irregular broad step lines produced by solution. 
The kinks in these broad step lines have again triangular or probably hexagonal 
contours which may be due to some local growth taking place at these step edges. 


§7. GROSS FEATURES 

The elementary spiral and other molecular growth features were only rarely 
found on the haematite crystals studied. The majority of the growth features 
observed were gross, consisting of triangular pyramids (fig. 10) with broken and 
irregular steps of heights of a few microns with a flat triangular top with rounded 
corners. ‘Three crystallographic directions inclined to each other at 120° 
appear as three directions of differential growth along the lines joining the three 
corners of the triangular feature to its centre. Similar observations have been 
reported previously (Maurice 1932). 


§8. COMPARISON WITH SILICON CARBIDE CRYSTALS 

Unlike the spirals on silicon carbide crystals, in haematite the point of 
emergence of the screw dislocation on the crystal surface, whence the spiral 
originates, is not marked by a hole, showing that in this case the core of the dis- 
location is not hollow. According to Frank (1951) a dislocation of Burgers’ 
vector exceeding about 10 is only in equilibrium with an empty tube at its core, 
which gives rise to a crater at the point of emergence of the dislocation on the 
crystal surface. ‘This conclusion will follow when the large strain energy near the 
dislocation is insufficient to overcome the surface energy which tends to close up 
thetube. However, for crystals where a complete dislocation is readily dissociated 
into a cluster of weak partial dislocations the core will not be hollow. Such is the 
case for the dislocations 47A high observed on paraffins by Dawson and Vand 
(1951). In haematite this does not seem to be a possible explanation. It is 
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likely that the magnitudes of the surface and strain energies depending upon the 
elastic constants of the crystal are such that the core closes up. 

Whereas in most cases of silicon carbide the spirals are more closely spaced 
near the centre than further out, in the spiral of fig. 1 the separation between the 
successive arms gradually decreases with increased distance fromthe centre. ‘This 
is attributed, not to the gradually decreasing supersaturation, but to the limitation 
placed by the triangular boundary enclosing the spiral. Further, in accordance 
with the theory, the step lines of the polygonal spirals on silicon carbide are observed 
to have a curvature near the centre, whereas in fig. 1 the step lines are observed to 
be quite straight and the corners quite sharp at the centre. The gradual rounding 
off of the corners and of the step lines takes place only away from the centre. 
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ABSTRACT. The concept of a kink-band is introduced. It is a thin plate of sheared 
material in a crystal, transverse to a slip direction, bounded by opposite ‘ tilt walls’ of 
dislocations. It is shown that if the angle of kinking shear y exceeds a critical value 
(estimated to be 34°) the stress concentration at the edges of the band is sufficient for the 
creation of new dislocation pairs or loops, and that if the applied stress exceeds a certain 
value, inversely proportional to the square root of the diameter of the kink-band, the kink- 
band will grow by this process. ‘The stress field of a kink-band is similar to that of a slip- 
band, a plate of deformation twin, or a crack under tangential stress. In suitable 
circumstances these may act as sources for each other. ‘The successive generation of 
kink-bands, inside each other, from the same source, will produce a kink of the kind 
experimentally observed. 
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§1. INTRODUCTION 


ROWAN (1942) has been foremost in maintaining that there is a mechanism 
Ca) of deformation of crystals, called by him ‘kinking’, which is quite 


distinct from the other mechanisms, slip and twinning. Speaking in © 


terms of dislocations, the essential difference is that slip requires the generation 


of a succession of many dislocations all in the same slip plane, while kinking 
requires the generation of dislocation pairs in many parallel slip planes, regularly | 


spaced a small number of atomic distances apart. In this paper we construct a 
somewhat idealized ‘kink-band’ in an otherwise perfect crystal subjected to a 
uniform applied stress, and show that the kink-band is able to act as a sufficiently 
effective concentrator of stress to extend itself by generating new dislocation 
pairs at its edges. 

The kink-band we consider is a region between two approximately plane 
parallel ‘walls’ of edge dislocations, of the kind considered by Burgers (1939, 
1940), Bragg (1940), van der Merwe (1950), and Read and Shockley (1949, 
1950). The dislocations in these two walls are of opposite sign, so that the walls 
are forced in opposite directions, away from each other, by the applied stress: 
but the walls are of finite extent, and their edges attract each other. 


§2. FIRST APPROXIMATION : PLANE WALLS 


It will give us an instructive first approximation if we first consider an artificial © 


constraint to be present, which keeps the walls strictly plane. The stress field 
of such a wall, of infinite extent, falls off with distance x from the wall as e 


p 


Fig. 1. 


where D is the spacing between dislocations. It is negligible at distances large 


compared with D. A finite wall of ‘height’ L measured normal to the dislocation — 


lines (infinite along their length) has a stress field like that of a large dislocation 
at distances large compared with L, diminishing inversely with distance; but 
at distances from either edge small compared with L and large compared with D 
the stress is independent of distance and depends only on the angular coordinate 
measured from the nearest edge. If (fig. 1) 0, and 6, are angular coordinates 
measured from either edge (0, =0=6, being the plane of the wall) the shear 


—a/D 
>a 
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stress on planes parallel to the glide planes of the dislocations, at any distance 
greater than D from the wall, is 
T oy = GA(sin 20, +s8in20,)/4a(1—oc)D, =. ae (1) 

where A is the strength of the dislocations, G the rigidity modulus and o Poisson’s 
ratio. This formula may be obtained by approximating the sum of stresses 
arising from individual dislocations by an integral. In consequence of this 
stress, a second wall of dislocations of the opposite sign, parallel and level with 
the first, at a distance w from it, small compared with L, is attracted with a force 


ROMmin Lieinll—o)D?-) 2) = Mak (2) 


per unit length measured along the dislocation lines. This result is likewise 
obtained by replacing the sum of the forces on individual dislocations by an 
integral. The maximum force is experienced by dislocations a distance 
approximately w from each edge, i.e. this attraction acts predominantly between 
the edges of the walls. Since it increases very nearly proportionately with w, 
the constant force tAL/D exerted on either wall by an applied shear stress 7 
will bring the walls into a position of equilibrium, defined by 


(w/L) In (L/w) =2(1—c)rD/GA ee) 


i.e. in order of magnitude, 


Ae SA GND Ike pet Se hae (4) 


The resultant stress, due to the superposition of the applied stress and the two 
stress fields of the form (1), which tend to cancel each other except on the mid-plane 
between the walls, has a maximum at points on this plane from which the two 
neighbouring edges of the walls subtend an angle of 90°, i.e. a distance w//2 
beyond the edges of the walls. This maximum stress is 


Peers GA eae o) ne wl Was Ome it. « (5) 
and hence, when the equilibrium condition (3) is satisfied, in order of magnitude 
Tmax ~7(1 ae L/w). aye ciel sie (6) 


Supposing that A/D (the angle of kinking) is of the order of magnitude 10-1, 
Tmax/G is of similar order of magnitude, and thus sufficient for the generation 
of new dislocation pairs; we must leave the question of whether this maximum 
stress occurs at the appropriate distance from the dislocations already existing 
till we have developed a better approximation, in §3 below. We infer that if the 
angle of kinking is sufficient, the kink-band can grow by making new dislocation 
pairs at its edges; it can of course also shrink by mutual annihilation of pairs 
at the edges. T'o find which will happen we must discover which process lowers 
the energy. 

The stress field of the pair of parallel tilt walls, observed at a distance large 
compared with D, is equal to that of two short walls bridging the gaps at the 
ends, as we may see most simply (fig. 2) by the consideration that the addition 
of two such walls of opposite sign would reduce the whole stress field to that of 
a slab of rotated crystal bounded by a rectangular parallelepiped of tilt walls, 
with a stress field effectively confined to a distance of the order of magnitude 
of Dfromthem. Seen at distances further away than a, this stress field is virtually 
that of two strong dislocations of opposite sign, their Burgers vectors being 
wA/D, directed in the plane of the main walls. ‘The energy of the system as 
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compared with that of the unkinked material under the same applied stress 
is approximately 
U =2G(wdA/D)? In (L/w)/477(1 —o) 
+2GA?(L/D)1n(D/1r1)/40r(1—0)—Lwra/D, wa ae (7) 
per unit length measured along the dislocation lines; 7, is a length, of the same 


order of magnitude as A, dependent on the ‘core energy’ of dislocations. These — 


three terms represent respectively the energies of macroscopic and microscopic 
stress fields of the dislocation array and the energy released by removing the 
previous state of strain from the kinked region. If, under constant 7, w/L retains 
its equilibrium value according to (3) this becomes 


GXLIn(D/X) (1-0)? L? 


"Ball =o)D' 26 ln (Loy? sei va as (8) 4 


which is a maximum, as a function of L, when 
L  G*aln (L/w) In (D/A) 
iment ee 
Hence the kink-band will grow if it has more than a critical depth, which is 
inversely proportional to the square cf the stress, and for A/D ~ 1071, A~ 10-8 cm 
and 7/G~ 10°, is of the order of 10-° cm. 
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Fig. 2. (a) may be regarded as the superposition of (b) and (c): (6) has no remote stress field 
(representing an-enclosed grain of crystal with different orientation): hence the remote stress 
field of (a) is the same as that of (c). 

This simple treatment has brought out the main qualitative features of the 
situation, which we shall confirm by a more accurate but less direct method. 
The degree of success of the simple treatment depends on the fact that the walls 
have considerable stiffness, and so do remain approximately plane. Nevertheless, 
because the attraction between them is predominant near their edges, we must 
expect these edges to bend in towards each other: we shall in fact conclude 
that in the equilibrium configuration they bend round to become the sides of 
a thin elliptic cylinder. 


§3. SECOND APPROXIMATION: THE ELLIPTICAL KINK 
Consider an infinite isotropic elastic body containing a cavity in the form 


of an infinite cylinder, whose generators are parallel to the z-axis of coordinates _ 


and whose cross section is the ellipse 


wlartyb=1, ba reeves (10) 
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The displacements and stresses under an applied shear stress 7,,=7 have been 
calculated by Starr (1928). We shall use his results, omitting terms which are 
of second or higher order in stress-dependence. 

Suppose the cavity had been filled with a jelly of elastic moduli negligible 
compared with those of the external material; the applied shear stress would 
have produced in this jelly a uniform shear 


ieee sie ae (11) 


Supposing the external body and the jelly both to be crystalline, with the same 
lattice structure in the unstrained condition, there is still atomic correspondence 
at the boundary of the ellipse in the strained condition. Now, suppose the jelly 
removed and replaced by unstrained crystalline material, identical with the 
external material. Since it is unstrained it produces no macroscopic stresses, 
and macroscopic elastic equilibrium still exists, but there is now atomic misfit 
at the interface. ‘This gives rise to local centres of strain, which may be identified 
as dislocations; the strains due to these are only appreciable at distances from 
the interface of the same order of magnitude as the distance between neighbouring 
dislocations. The dislocation density depends on the orientation of the plug of 
unstrained crystal inserted into the cavity, and is a minimum when this is chosen 
so as to give continuity of lattice planes across the ‘thickness’ of the elliptical 
plate. ‘The dislocations present then all have their Burgers vectors directed 
along the x-axis, and their distance apart in either ‘wall’, measured along the 


y-axis, 1S D=Aly=Gadjrb(1—c), nena (12) 


where A is the ‘strength’ of the dislocations. 

Our calculation is analogous to that of Leibfried (1951) for finding the 
equilibrium distribution of various dislocation systems all in one slip plane, 
under various stress fields. Several problems of this kind were treated in terms 
of discrete dislocations by Eshelby, Frank and Nabarro (1951). Leibfried’s 
treatment of the same problems by continuum theory leads to results which are 
asymptotically identical for the case of large numbers of dislocations, and very 
similar, displaying all the principal qualitative features of the equilibrium 
distribution, even for small numbers of dislocations. We infer by analogy that 
the same applies to our present problem, and that we have thus discovered, to 
a good approximation, an equilibrium configuration for the dislocations bounding 
our idealized kink-band (fig. 3). They lie on the surface of an elliptic cylinder 
of cross section defined by (10) and have the spacing D defined by (12). 
Rearranging our equations, we may say that for a constant dislocation spacing D 
(or constant kinking shear y=A/D) and a constant depth of kink-band 26 the 
thickness of the kink-band 2a varies with the applied stress 7 according to the 


Bayon 2n=2br(1—eV/Gy HOF grin AE (13) 
in agreement with (4). 
The radius of curvature at the edge of the kink-band is 


pa Dba lmoyYG yee ackeee (14) 


The stress field around the kink-band, at points not too close to individual 
dislocations, will be that found by Starr. In particular the shear stress 7,,, has 
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a maximum at a point opposite each end of the ellipse and distant p from it and 


is there equal to ro ueleis a aiinla 
which is comparable with (6). 


If y is constant, then by eqn. (11), b7/a is constant, and hence ty, is 


independent of 7: 


Tmax = (2/3 »/3)Gy/(1 —¢) oat ats (16) 
(ct, eqnin(S))a lt Tree Tel). 2g ye Cee ee eee (17) 


the yield stress of the perfect crystal it should generate a new dislocation pair 
and so deepen the kink-band. A kink-band can therefore extend by the proposed 


mechanism if it has not less than a certain characteristic kinking shear y, 
determined by natural constants of the material: 


¥ 2¥_= 3/371 —«)/2G. 


If we take o=4 we have the simple formula 


Vs Brel Gay iO her Oe ee (19) 


From Mackenzie’s estimate (1949) appropriate to inert gas crystals and, for 
want of a better model, to metals, +, ~ G/30, giving y,.~34°. 
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Fig. 3. Diagram shows the atomic configuration and the dislocations 1n a kink-band. 


If the angle of kinking exceeds the critical value the kink-band can grow by 
making new dislocation pairs at its edge; it can in any case also shrink by © 
annihilation of dislocation pairs. Which process will happen is decided by which 
lowers the energy. 

The decrease of energy in the kinked configuration, relative to the unkinked 
configuration under the same applied shear stress, is obtained by taking Starr’s 
(1928) value for the energy decrease (the formula for [W] on p. 500 with sign 
reversed, neglecting the small quantity A(&)), and adding the ‘microscopic’ 
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| strain energy of dislocations in the surface. This gives, for unit length along the 
| dislocation lines, 


= —n(1—0)1°b?/2G + GbdyIn(A/yr,)/m(1—0), cee ee (20) 


where 7, is a length of the same order of magnitude as A, dependent on the 
imperfectly known ‘core energy’ of a dislocation. U is a maximum, with 


respect to 6, when b/A=G?yIn(A/yr;)/r2m2(1—o)® nas (21) 
aC eee ee eee (22) 


Inserting numerical values, the smallest kink-band which will grow under a 
stress of 10-°G is about 10-3 cm deep and 10° cm thick. 

We now bring in the requirement that the kinking should continue in the 
Same manner, with a constant kinking angle y. The second term in (20), 
expressing the local energy of dislocations is represented as varying continuously 
with 4: actually it will vary periodically, as new dislocation pairs are formed. 
We suppose that the configuration pauses at positions where the last pair of 
dislocations are a distance D/2 from the end of the ellipse, this being the 
configuration in which the strain due to the dislocation array corresponds most 
nearly to continuous kinking at the elliptical boundary. ‘Then for uniform 
continuation we require the next pair of dislocations to be created a distance 
D/2 beyond the edge of the ellipse. If we equate this to p, which gives the 
position of maximum stress according to Starr, we have, from (14) 


ID=AQy=p=Cl =o)br/Gs?, aces (23) 
whence G(T 6) NyAfzpee te (24) 


which is (1 —«)"1, i.e. about 14, times the stress required according to (22). 

We commit a logical error here by employing Starr’s formula too close to 
the dislocation system, at a distance of the order of D instead of a distance much 
greater than D. However, the same conclusion is approximately confirmed by 
considering discrete dislocations. The position of maximum stress should not 
be far from the position of maximum stress from the nearest pair of dislocations, 
which is on the mid-plane between them at a point where they subtend an angle 
of 60°. If they lie on the ellipse, a distance H in the y direction from its end, 
their separation is H14#b'?723/2(1 —o)/yG, and if H is D/2, the point of maximum 
stress due to them is a further D/2 beyond the end of the ellipse when 


er aa MU ALND Ug eae aly (25) 
We reach much the same conclusion if we require that this point of maximum 


stress due to the last pair of dislocations should lie on the tip of the ellipse when 
they are D above it: this gives 


HORT sey Aeneas Ae ne. (26) 


Hence we conclude that there is a shear stress, dependent on the size of the 
kink-band as in (22), and not very much larger than is given by (22), at which 
the kink-band propagates with constant kinking angle. 

It is necessary to enquire whether the propagation at constant angle is stable. 
It appears by the argument above that an insufficient stress can allow the kink-band 
to propagate a certain distance but with increasing kink angle, while an excessive 
stress could produce dislocations too far away, causing the kink angle to diminish 
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until it is too small for propagation to continue. We may, however, use the 
following argument which implies that there will in general be a stable condition, — 

The condition 7 >7, defines a region S round the point of maximum stress 
within which the new dislocation pair will be created. In this region there will | 
be a continual production of dislocation loops due to thermal fluctuations, 
but in most cases these will disappear in succeeding fluctuations. For a loop to 
grow under the applied stress it must be carried by the fluctuations past some 
critical size. Accordingly we introduce a length /, characteristic of the material | 
and of the order of a few atomic spacings, and require that S should have linear 
dimensions not less than / before a new pair of dislocations may be formed. 

As a first approximation we may assume that the shape of the kink-band 
near its end does not change as the kink-band widens under the applied stress. 
Then the size of the region S and its distance from the end of the band will 
increase in constant ratio k to one another. We have stipulated that S shall | 
have dimensions / when the new dislocation pair is formed; it follows that this 
pair will be formed at a distance //k=h, say, from the end of the kink-band. | 

We have now to consider how h will depend on the spacing D of the 
dislocations in the kink-band. 


D— 


Fig. 4. Position of formation of new dislocation pair, against spacing in kink-band. 


If D>D)=4/y.. no new production of dislocations can take place. If D is 
slightly less than Dy the region S will be relatively small, that is k is small or 
h large. As D decreases the relative size of S steadily increases and h decreases. 
If h became zero for any value of D it would mean that the new dislocation pair 
was formed on the same glide plane as, and between, the end pair of dislocations 
of the kink-band. But the stress field of the latter is such as to oppose the 
formation of new dislocations, and since these dislocations are the most favourably 
situated with regard both to distance and to orientation factor we can expect 
them to contribute an appreciable part of the total stress field. Consequently — 
we do not expect that h can decrease as much as to zero for any D. ‘These 
considerations lead us to expect a curve of h against D of the form shown in > 
fig. 4. ‘The line h=D either may not meet this curve, or it may intersect it in two | 
points as shown in the figure. In the former case / is always greater than D, 
so that as the kink-band develops the spacing becomes larger until it eventually — 
exceeds Dy when no further dislocations can be produced. In the second case 
where h=D meets the curve the points of intersection give two values of D, 
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} D, and D, (D,>D,) which can be continued uniformly. It is easily seen that 
D, is stable and D, unstable. Hence if the spacing in the kink-band is initially 
less than D, it will tend to the value D,, while if the spacing is greater than D, it 
will increase until it exceeds D) when further growth of the band is prevented. 

This argument is based on continuum theory. When attention is paid to 
the existence of a crystal lattice, we must notice that there are certain preferred 
angles of tilt, those for which the mean number of lattice planes between 
successive dislocations is an integer, which give cusped minima of surface energy 
with respect to angle of tilt (see Read and Shockley 1950). It is possible that 
stable kinking could set in on one of these preferred angles when continuum 
theory would not predict stability. 


§4. DISCUSSION 

We have dealt with a kink-band in the form of an elliptic cylinder for reasons 
of simplicity, in order to have a two-dimensional elastic problem, with some 
part of the calculation done for us in advance. In physical reality we are more 
likely to have kink-bands in the form of ellipsoids, roughly penny-shaped. We 
suppose that this case would give us more work without teaching us anything 
new in principle, or even making much quantitative change in the results. The 
remote stress field of the penny-shaped kink-band will be virtually that of a 
closed loop of dislocation line, roughly following its edge, but with somewhat 
smaller diameter, with its Burgers vector proportional to the thickness of the 
band. 

Bounded slip-bands, kink-bands, plates of deformation, twin or martensitic 
transformation product, and thin cavities, cracks or notches under tangential 
stress, all produce stress fields of the same kind: in various circumstances any 
one of these may serve as the source of any other. As an example, in a substance 
such as rock-salt with slip planes, and slip directions at right angles to each other, 
a slip-band held up at any point might continue as a kink-band employing slip 
on the perpendicular (110) plane. The kink-bands observed by Brilliantow and 
Obreimow (1934) in NaCl may have such an origin. In the case of zinc a suitable 
initiator of a kink-band would be a small included twin of the main crystal, 
mechanical ‘untwinning’ of which would occur under a compressive stress 
along the base-plane. 

Experimentally observed kinks are usually much thicker than the thin 
kink-bands of our theory. However, it is to be expected that: thin kink-bands 
will build up into a broad kink. When a kink-band extends as far as a free surface, 
the attraction between the edges of the tilt walls disappears, and they should 
become parallel planes. A continuing stress should then force the walls apart. 
Now, suppose the kink-band started from some permanent source, such as a 
cavity or plastic inclusion. ‘The whole force acting on the walls is available to 
detach them from this source and, as soon as they separate from it, it should be 
able to initiate a second kink-band within the first. The sequel, leading to the 
‘stove-pipe’ kink which 1s typically observed, is shown in fig. 5. 

The last figures in the sequence show another phenomenon: the union of 
two tilt walls into one. This can occur if the component dislocations are not 
on the same slip planes, and in general lowers the total boundary energy. 

It should be noted that the slip planes on either side of a wall of a kink-band 
are symmetrical with respect to it. The slip planes inside and outside the band 
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(d) 

Fig. 5. (a2) Specimen with included twin now ‘ untwinned’ (or other stress raiser, e.g. cavity) 
showing stress concentration regions at its ends, and horizontal slip planes; (6) a kink-band 

has developed; (c) it has reached the surface; (d) its walls have moved apart, and a new 

kink band has commenced; (e), (f), (g), (4) the process continues: at (g) the outer walls 


have partly united. 


(c) 


Fig. 6. (a) Incipient deformation band, due to slip held up at an obstacle: O, O, regions of 
magnified and reversed stress; (b) production of transverse kink-band on reversal of applied 
stress (assuming the initial stress does not reverse); (c) the growth of a deformation band. 
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make an angle y with each other, so that the kink-band is tilted at an angle $ y 
from the normal to the slip direction in the external material. 

Figure 6(a) shows an incipient deformation band: a region of S-shaped 
lattice curvature where slip on slip planes 5S S is held up against some unspecified 
obstacle (any continuous slip through the specimen which may have by-passed 
the obstacle is irrelevant, and is omitted from the diagram). ‘This accumulation 
of dislocations of opposite sign is somewhat similar to our kink-band, and 
produces regions of magnified (and reversed) sHear stress at the points O. This. 
may be sufficient to produce new dislocation pairs, and so the commencement 
of a kink-band, transverse to the slip plane, but this kink-band will not grow 
indefinitely because the slip in it is in the reverse direction to that directly 
favoured by the applied stress. It could, hypothetically, develop into a simple 
kink-band across the specimen on reversal of the stress, as shown in fig. 6 (0); 
but under continuing stress of the same sign only a limited amount, if any, of 
pair creation will occur. In any case, it does not destroy the region of reversed 
stress. Hence dislocations on subsequent slip bands will be held up at the edge 
of the region of reversed stress. ‘Thus it will propagate across the whole specimen 
as indicated in fig. 6(c) and observed by Cahn (1951). The stress field of a 
developing deformation band is thus very similar to that of a kink-band. The 
differences are that the dislocations come from outside instead of from inside, 
so that the kinking is in the reverse direction, and are in groups side by side on 
the same slip planes, so that the lattice planes are curved, instead of sharply bent,. 
at the walls. ‘This is a situation of higher energy, and on heating ‘ polygonizes’ 
to something more nearly resembling the simple kink. Figure 6 seems to suggest 
that the initial obstacle must extend right through the specimen, normal to 
the plane of the paper, but this is unnecessary. The lattice curvature will 
spontaneously spread out in this direction from any small obstacle. In terms. 
of dislocations this is a consequence of the repulsion between the similar screw 
dislocations which unite pairs of opposite edge dislocations on either side of 
the band, but its necessity is seen most simply by trying to force an S-shaped 
bend locally and not elsewhere in a sheet of paper. 

There is more that could be said, for example with reference to twinning 
and martensitic transformations, and to creep, and possibly fracture, of 
polycrystalline metals, but we think we have said enough to indicate the 
considerable range of application of the theory of kinking, beyond the pure 
oe of kinking itself. 
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LETTERS 2O, THE EDITOR 
Supplementary Note on Ray Tracing 


It seems desirable to record two small additions to the original descriptions of the — 
methods of tracing rays through a symmetrical optical system (Smith 1920, 1945). Accord- 
ing to the first of these methods the difference of path for a ray lying in an axial plane, which 
is refracted at a spherical surface, is represented by a fraction, the denominator being 


D=4 cos 44 cos $¢ cos $¢’ cos $y’. 


As the cosines of these half angles have not been found, this expression is far from convenient. 
The cosines of the complete angles, as well as the sines, are however already known with 
ample accuracy. In general the most convenient equation to use is D=S-+P/S, where S _ 
is the sum of the four cosines and P is the product of the four sines. ‘The second term 
in this expression for D is often negligible, frequently makes a contribution to the fourth 
decimal place only, and only occasionally amounts to a third decimal place correction to S. 
In tracing rays by the second method (which is available for tracing any ray through 
any surface of revolution) I learn that W. Weinstein (1952) has found cases where the 
published procedure fails to converge. The trouble arises in finding a coordinate of the | 
point of intersection of a ray and the surface when the equations of both are known. The 
equation of the surface is taken in the form s=f(z) where s is the subnormal and 2 is the 
coordinate measured along the axis of symmetry. The calculation begins by guessing a 
value of x, say 2%. Then the point xo, yo, % on the ray will also lie on the surface if 


Xo? +o? =2 J sdz. 
If this condition is not satisfied a new value of z is found, viz. 21, given by 
25921= Xo +o +2 fz ds, 


where sy is the subnormal for points in the plane z=2, and the limits of integration are 
sufficiently obvious. The same procedure is followed with 2, as the assumed value of z 
instead of 2), and this yields a further value of z, say g,. In nearly all cases 2, 21, 20,..-- 
is a rapidly converging series leading to the value of z required. But in some exceptional 
‘cases, where the ray makes a large angle with the axis of symmetry, it may diverge. ‘The 
question to be considered is how the calculation should then proceed. The answer is 
simple: 2, is not the best approximation that can be taken for the continuation of the 
computation, though in almost all cases it has been found so good that no improvement has 
been required. The value now suggested for further approximations will be denoted by 2’. 
Take 2, as a reference point for the measurement of the differences between successive 
approximations, and write 


Z%=2,+d, and Zg=21 +d. 


If the values found according to the original scheme are rapidly convergent, d, will be much 
smaller than d,. But in any event a better approximation than 2, is 2’ where 


d, -d, 
‘The additional labour involved in finding 2’ is clearly trivial, for optical workers regularly 
use tables of reciprocals suitable for this application. Alternatively z’ is given directly in 
terms of 2, 21, and 2, by the equation 2’=(z92,—2,”)/(%)—221+22). It is thought that 
when this modified routine is adopted the point in which the ray meets the surface can 
always be found rapidly with whatever accuracy is required. 

A somewhat extreme example illustrates this, eight decimal places being retained to show 
the efficiency of the modified procedure clearly. 

The ray x=0, 0:28y+0-962=1:104 meets the paraboloid x?+y?=1-8z in the point | 
y=1-2, z=0°8. Suppose the approximate method is employed with the original guess 
%)=0°81. The next two approximations are 2,;=0-75493878, 2,=1-01925467, showing 
that the original procedure leads to a strongly divergent series. These values of 2, 21, 25 


2’=2,+d'’ and 


once, yo 
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give zo =0-80050734, and treating this as a starting value the regular routine yields 
2, =0-79768241, z,’=0-81062976. Applying the new formula again to 29’, 21’, 22’ gives 
the improved value z)”=0-80000089. The series of values of z found by using the modified 
routine clearly converges rapidly to the correct limit. 


Roselyn, Holton, T. SMITH. 


Wincanton, Somerset. 
30th June 1952. 


SMITH, I. 1920, Proc. Phys...S0c.,.32, 2525 1945, [bid., 57,286. 
WEINSTEIN, W., 1952, Proc. Phys. Soc. B, 65, 731. 


Positive and Negative Joshi Effect in a.c. 
‘Silent ’ Discharges in Iodine Vapour 


It is well known that the Joshi effect is associated with current pulses which can be 
observed when the potential difference across a suitable high resistance, placed in series 
with the a.c. discharge circuit, is delineated on the fluorescent screen of an oscillograph. 
The object of the present communication is to show some typical oscillograms obtained 
with an iodine vapour ozonizer tube for different applied a.c. voltages which indicate two 
distinct sets of current pulses, and to draw certain definite conclusions regarding their 
nature and origin. Some experimental results suggesting different origins for the two 
distinct sets of pulses are also described. 

Figure 1 shows a few typical oscillograms obtained with an iodine vapour ozonizer, 
alternately exposed to and shut out from light, for a range of 50 c/s voltages. These oscillo- 
grams showed the following features with successively increasing applied voltages: (i) Starting 
from a low applied voltage, and critically adjusting it so as to secure just the sinusoidal 
50 c/s current trace in the dark, it was found that distinct pulses appeared on irradiation, 
each as a line on the 50 c/s trace near each peak. The critical voltage in the particular case 
was 160 volts. (ii) At slightly higher voltages these pulses appeared, even in the dark, with 
an increase in number and height and persisted on irradiation. (iii) At a slightly higher 
critical voltage (200 volts) another set of longer pulses appeared in the dark, one, and 
sometimes more, near each peak of the 50 c/s current, together with the previous set of 
pulses. The two sets of pulses showed a distinct difference in height, position, number 
and distribution. On irradiation, the longer pulses were either quenched altogether or 
reduced considerably in height, and the shorter pulses were found to persist with an increase 
in number and sometimes in height. (iv) Both sets of pulses observed in the oscillograms 
taken in the dark increased in number and height with increasing values of the applied 


e voltage. (v) The reduction in height of the longer pulses on irradiation was found to be 


gradually smaller as the applied voltage was increased. (vi) At high voltages the two sets 
of pulses could not be distinguished from each other, and the pulse patterns remained 
practically unchanged on irradiation. 

We are inclined to think that the longer pulses in our oscillograms should be considered 
as due to discharges taking place between the negative surface charge on the glass surface 
of either electrode and the positive ions in a manner first suggested by Klemenc, Hinten- 
berger and Héfer (1937), and later reported in a slightly modified form by Deb and Ghosh 
(1946, 1948) at the suggestion of Mitra. This view has been the basis of an explanation 
of the negative Joshi effect. 

The other set of comparatively short pulses observed in the oscillograms should be 
identified with the electron avalanches formed by the Townsend collision process. With 
the voltage critically adjusted to obtain just the 50 c/s trace in the dark, the pulses which 
flash out on irradiation are considered to be due to the formation of avalanches by the 
photo-electrons released by the layer of negative surface charge on irradiation. With 
slightly higher voltages such Townsend pulses will be formed, even in the dark, by the 
electrons in their flight towards the anode in the normal way.* In the case of an absorbing 
gas or vapour like iodine it is also possible that the absorbed light energy may provide 
energy for accelerating the electrons sufficiently to produce Townsend avalanches. 


* This view regarding the origin of the positive Joshi effect will be given in full elsewhere. 
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The two distinct sets of pulses—the ‘discharge’ pulses and the Townsend pulses— 
were also observed with an iodine vapour discharge tube fitted with adjustable external 
‘sleeve’ electrodes. On examining them with a high sweep-frequency of the linear 
time-base of the oscillograph, each pulse of either set was found to have either a continuous 
or an oscillatory form of discharge according to the internal resistance value of the discharge 
tube. The continuous form of the Townsend discharge pulses which appeared on 
irradiation for a certain distance between the two electrodes when the applied voltage 
was near the ‘threshold’ value is shown in fig. 2. It can be seen that in each discharge 
the current at first rises, attains a peak value and then falls rapidly with time. With an 
increased distance between the electrodes these pulses were found to change over to the 
damped oscillatory form. At a certain higher voltage the ‘discharge’ pulses appeared 
in the dark as well as the Townsend pulses. By adjusting the distance between the electrodes, 
and by using a higher sweep-frequency, it was possible to observe the two sets of pulses, 


Fig. 1. Pulses in iodine vapour ozonizer 
at different a.c. voltages (50 c/s). 
The oscillograms with  light- 
shutter closed are on the left and 
those with light-shutter open are 
on the right side. 


Fig. 2. Townsend pulses in 
iodine vapour ‘sleeve’-dis- 
charge tube observed on 
irradiation at acomparatively 
low applied a.c. voltage with 
a sweep-frequency corres- 

. : ponding to 1500 micro- 

Big. 1e seconds. 


one as continuous and the other as oscillatory, or both as continuous or oscillatory 
discharges having distinctly different damping coefficients. The ‘discharge’ pulses 
could be easily identified as they disappeared on irradiation at some suitable applied 
voltage. The different origin of the two sets of pulses is strongly suggested by these 
observations. The duration of both sets of pulses was found to be of the order of 
10~* to 10-4 second. 

It should be noted that the photo-suppression or the photo-reduction of what has 
been called the ‘discharge’ pulses is considered here as the main cause of the negative 
Joshi effect. It is, however, worthy of notice that the theory put forward by Joshi (1939, 
1945, 1946) can explain in a general way the reduction of Townsend pulses on irradiation. 
The photo-dissociation theory of Harries and von Engel (1951) can also explain the same 
in the case of chlorine. It seems to us that these physical processes may be operative in 
addition to the photo-suppression of the so-called: ‘ discharge’ pulses under suitable 
experimental conditions. 


Wireless Laboratory, S. R. Kuastecir. 
Physics Department, College of Science, P. S. V. SeTry. 
Benares Hindu University. 
16th July 1952. 
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The Epitaxial Growth of Germanium on Rocksalt 


A study has been made of the conditions necessary for the epitaxial growth of 
germanium on the (100) cleavage face of rocksalt with a view to using the method for 
preparing thin, single-crystal films of this material. Kénig (1944) has examined films of 
germanium deposited on KBr at room temperature and finds them to be amorphous. 
Films which had been heated (after removal from the substrate) to 460°c—500°c were 
found to be polycrystalline, with large crystallites in random orientation. Substantially 
similar behaviour has been observed in these experiments for films deposited on rocksalt 
maintained at temperatures up to about 480°c. 

The germanium was evaporated from a tantalum boat; to avoid the possibility of 
atmospheric contamination, the films were prepared in the specimen chamber of the 
electron diffraction camera at a pressure not exceeding 510-5 mm Hg. No trace of 
oxide could be detected in films so produced. Experiments were carried out at different 
rates of deposition of the germanium, ranging from 5 A/sec to 100 A/sec. The variation of 
beam intensity was effected by adjusting the evaporation temperature over the range 
1100°c to 1300°c. Although this procedure alters the mean velocity of arrival of the 
atoms as well as the beam intensity, the small change in velocity over this range of 
temperature is unlikely to influence the deposition process. “The presence of the electron 
beam during the condensation of the germanium was found to inhibit the deposition of 
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the film, the thickness of the film on the part of the target traversed by the electron beam 
being only a few per cent of that on the neighbouring parts of the surface (fig. 1). "The 
structure of the film is, however, unaffected by the presence of the beam. 

The first signs of an orienting influence of the substrate on the germanium occur at 
a temperature of 490°c. Sharp spots appear, superposed on continuous rings, indicating 
the presence of some crystalline regions with (100) and (111) faces parallel to the substrate. 
For (100) orientation the cube edges of the germanium are parallel to those of the rocksalt. 
For (111) orientation both arrangements in which the cube face diagonals coincide are ~ 
found. ‘These regions probably develop from cleavage steps on the rocksalt surface, at 
which the substrate influence is greater than on the plane facets. (Examination of cleaved 
rocksalt faces by multiple-beam Fizeau fringes shows that such surfaces are crossed by 
many such steps, a few tens of atoms in height and running mainly in [010] and [011] 
directions.) At this temperature deposition can be effected at low and high rates of 
evaporation. The films formed at the low rate show the reflectivity characteristic of the 
bulk material; those formed at the high rate are grey and lustreless and are seen to consist 
of large irregular aggregates. On raising the temperature to 510°c deposition is found to 
occur only at the high rate of evaporation, due to the appreciable evaporation of the rocksalt 
at this temperature. Fairly complete (100) orientation of the film occurs (fig. 2), only a 
trace of randomly oriented material being present. The degree of orientation is found to 
increase as the temperature is raised to 530°c. Above this temperature the rate of 
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evaporation of the rocksalt is too high to allow formation of the germanium film at the 
highest rates of evaporation used (0-01 g/sec into a solid angle 277; source to target distance 
6 cm). The films examined ranged in thickness up to about 3000 A; no dependence of 
structure on thickness was observed. 

It was found to be possible to maintain a (100) rocksalt face supporting an epitaxially 
grown silver film at a temperature of 600°c; germanium films were deposited under 
this condition. It was expected that the increased surface mobility of the germanium 
atoms at the higher temperature would result in the formation of more perfect mono- 
crystalline layers. This was found not to be the case (the films being similar to those grown 
on rocksalt at 510°c), probably due to the different nature of the substrate forces present. 


University of Reading. ; L. E. COo.uins. 
18th July 1952. O. S. HEAVENS. 


Konic, H., 1944, Reichsber. Phys., 1, 4. 


REVIEWS OF BOOKS 


Glass and W. E. S. Turner, 1915-51: Articles by Twenty Contributors, edited by 
F. J. Goopinc and E. Mricu. Pp. 144. (Sheffield: The Society of Glass 
Technology, 1951.) 20s. 


-Emeritus Professor W. E. S. Turner, F.R.S., celebrated his 70th birthday on 22nd 
September 1951. He was a founder of the Society of Glass Technology in 1916 and was its 
Secretary from that date until 1945. He was editor of the Journal of the Society from its 
foundation until 1951. In recognition of his distinguished career, twenty ex-students 
contributed the chapters which form the volume under review. This volume constitutes a 
record, even if somewhat superficial in parts, of the expansion of glass technology in this 
country during the period of Professor Turner’s association with the industry. It should 
not be inferred from this that all the developments in the industry during that time have 
arisen from that association, but it is true to say that he had more influence on it that any 
other single person. The book has internal evidence of its variety of authors, and some of the 
contributions are more professionally written than others. This is probably unavoidable in 
a publication of this kind, and the two main editors obviously took the view that it was better 
to do the minimum of editing even though this would inevitably mean variations in style. 
This decision was emphasized by the urgency with which the volume had to be produced, and, 
as stated in the text, this is the reason why all references and diagrams have been omitted. 

The range of the volume, and by implication the range of Professor Turner’s activities, is 
considerable, covering not only glass technology but also the ‘encouragement of art and 
design’ and the ‘ evaluation of refractories for glass making’. ‘The volume was obviously 
not intended as a textbook for experts. It is written in a style which makes it attractive to 
those who have a general interest in the industry or in the application of science to industry, 
and in no place does it become so highly technical as to be likely to cause difficulty to the 
more general reader. ‘There are one or two points to which the attention of such a reader 
should be drawn, as the statements made may give rise to wrong interpretations. Page 38, 
for example, states that barium is used as carbonate or oxide; it is more usual to use the 
carbonate or nitrate, and the use of oxide must be extremely rare. A further statement on 
this page that barium flints are no longer manufactured commercially is also incorrect. 
On p. 58 there is a sentence which is obscure, where the statement is made that no mechanical 
gas-producer plant has been developed for the firing of pot furnaces. It seems to be a 
matter of indifference whether the gas is used for firing pot furnaces, tank furnaces, or any 
other type, as the point at issue is the production of a mechanical gas-producer plant, and 
these are well known. On p. 119, table 22, the figures for coefficient of expansion for glasses 
6 and 7, namely Hysil and Monax, have apparently become transposed, the lower figure 
referring to No. 6 and the higher to No. 7. These however, are minor criticisms in a publi- 
cation which, apart from its intrinsic interest, is a compliment to a great man, is interestingly 
and attractively written, and covers a very wide subject adequately ina relatively restricted 
number of pages. W. M. HAMPTON. 
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Coloured Glasses, by W. A. WEYL. Pp. xvi+541. (Sheffield: The Society of 
Glass Technology, 1951). 35s. 


This book is the fifth monograph on glass technology published by the Society of Glass 
‘Technology under the general editorship of Professor W. E. S. Turner. The range covered 
by the book can best be indicated by stating that it is divided into five parts under the 
following general designations : 

Constitution of’coloured glasses. 

The colours of glasses produced by various colouring ions. 

The colours of glasses produced by non-metallic elements. 

The colours of glasses produced by metal atoms. 

Fluorescence, thermoluminescence and solarization. 
In view of the range of subjects dealt with it is not to be expected that practical information 
immediately useful to the glass maker could be included, and it is clear that the author 
has written with an eye on the scientific as distinct from the production staff of the glass 
industry. On the other hand the approach is so specialized that it is unlikely the book will 
appeal to a much wider public than is provided by the glass industry, although the earlier 
portions, dealing with the constitution, will obviously be of value to physicists and others 
working in related fields. The book is very fully referenced and has two good indexes. It 
is well produced, with clear diagrams, but is, in the nature of things, somewhat difficult to 
read. ‘The main criticism which could be levelled against it is that while the author has 
covered practically the whole of the literature on the subject up to the year 1949, he has done 
so in a rather uncritical manner. Some of the papers which are reported are of a much 
higher standard than others and there is no indication or comment by the author as to whether 
the conclusions have been modified or not by later investigators. As long as the reader will 
remember that he must make his own critical appreciation of the value, then this series of 
extended abstracts will be of considerable use. The author states that the production of the 
volume covered some 12 years, and the amount of information that is collected into one 
volume could, one supposes, scarcely have been produced in less time. 

W. M. HAMPTON. 


Radio Astronomy, by BERNARD LovgLt and J. A. Ctecc. Pp. 238. (London: 
Chapman and Hall, 1952.) 16s. 


The rapid growth of radio astronomy since the war has made it difficult, even for those 
working in related subjects, to keep in touch with the various investigations which have 
developed from the original work of Jansky and others in the 1930’s. This difficulty is 
largely due to the diversity of the techniques and experience which are required in applying 
radio methods to astronomical investigation. Professor Lovell and Dr. Clegg have succeeded 
in giving a clear and authoritative account not only of the fundamental principles and 
possibilities of this new subject, but also of the experimental methods and results of many 
of the investigations which have already taken place. 

The publication of this book is particularly opportune in that the first exploratory stage 
of the subject is perhaps now over; many preliminary studies have been made and the 
directions of future methods of investigation are beginning to become clear. One line of 
approach is illustrated by the plans, recently announced, for a large rotatable radio telescope 
to be built in this country under Professor Lovell’s direction. The developments which 
have led up to these plans, and the experiments which should be possible with such an 
instrument, will be of interest to a wide group of readers. 

No specialized knowledge either of astronomy or of radio techniques is assumed by the 
authors, and excellent introductory chapters are given on the fundamentals of both subjects. 
A detailed account is given, in the following six chapters, of the applications of pulse radar 
methods to the observation of meteor trails. It is shown clearly how vital has been the 
contribution of these methods to meteor astronomy, not only in permitting daytime measure- 
ments but in providing more accurate data than has been available from visual observations. 
In particular the radar measurements of meteor velocities have succeeded in settling the 
important question as to whether a significant fraction of meteors is of interstellar origin. 
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The following chapters deal with the observations and interpretation of the radio emission 
from the sun and the galaxy; the possibilities of investigating the structure of the solar corona 
by radio methods, the effects of sunspots and of solar eruptions and what is perhaps the 
most intriguing of the present problems, the nature of the ‘ radio stars’. Short accounts 
are given of the ‘ twinkling ’ of radio stars due to irregularities of the ionization in the upper 
regions of the ionosphere, radar observations of the ionization associated with aurorae, 
and of the methods of investigating the moon by radar and radio methods. In a final 
chapter the possibility of measuring the distances of planets, and hence of determining the 
solar parallax, by radar methods is discussed. 

Throughout the book an excellent balance has been maintained between the descriptive 
and the mathematical treatment. It should prove thoroughly readable for the general 
reader with some scientific background, yet at the same time it contains sufficient funda- 
mental detail to satisfy those who wish to enquire more deeply into the subject. A good 
selection of suggestions for further reading is also given at the end of each chapter. M. R. 


Electrical Phenomena at Interfaces, edited by J. A. V. BuTLer. Pp. vii+309. 
(London: Methuen, 1951.) 32s. 6d. 


The volume under review is a successor to Electrocapillarity by J. A. V. Butler, wie 
appeared in 1940, and is very largely a new work in which the editor and six others have 
collaborated to produce a volume of very considerable interest. Electrical phenomena at 
interfaces ‘pop up’ with remarkable persistence in various fields of physics, chemistry, 
and biology, and it is valuable to have so many aspects of the phenomena brought together 
in one volume. The wide range of subject matter makes the book rather difficult to review; 
but careful editing and choice of topic have overcome most of the difficulties which 
inevitably arise from a collaborative effort of this nature. 

The book opens with two chapters by J. A. V. Butler on electrical potential differences 
and their origin, and the electrical double-layer and electrocapillarity. There follows a 
discussion by J. M. Creeth of electrokinetic effects, and an account by P. A. Charlwood 
of electrophoresis of proteins. ‘The electrokinetic behaviour of ionic crystals is discussed 
in a short chapter by 'T. R. Bolom, and the stability of colloidal solutions is analysed by 
M. B. M’Ewen, particularly in relation to the Verwey—Overbeek theories. Three chapters 
more in the domain of conventional physical chemistry follow : overpotential (J. O’M. 
Bockris), concentration polarization and the deposition of metals (J. A. V. Butler), and the 
behaviour of oxygen and hydrogen at electrodes (also by the editor). The last two chapters, 
by W. F. Floyd, bring out the importance of electrical phenomena at interfaces in biological 
systems. ‘The first of these chapters discusses membrane potentials and the electrical 
properties of living cells, and finally there is an account of the electrical properties of 
nerve and muscle. 

The book is well written throughout and there are many references to original papers. 
Altogether it should be of value to research workers and to those whose task is to teach 
particular features of this growing field of knowledge. J. T. RANDALL. 


Polarisation de la matiére. Pp. 171. 1st edn. (Paris: Centre National de la 
Recherche Scientifique, 1949.) 1800 fr. 


This is a collection of papers on dielectric and magnetic properties of matter, read at a 
conference organized by the French Centre National de la Recherche Scientifique. The 
main emphasis is on electric polarization, and the collection forms a useful review of modern 
theory and experimental methods on electric dipole moments. Ay HGe 


Prisms and Lens Making, by F. 'Twyman. Pp. viiit629. 2nd edn. 
(London: Hilger and Watts, 1952.) 58s. 


This second, enlarged edition of Prism and Lens Making is four times the size of the 
first and is a far more important contribution to optical literature. Mr. Twyman mentions 
in the preface that in writing it he has been assisted by contributions from other specialists, 
and the text contains many long extracts from original papers. ‘These add interest and 
vividness without destroying the unity of the book or submerging the author’s own agreeable 
style. 
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Chapter 1 is an historical introduction with many interesting quotations from original 
sources. It does not, however, include Leonardo da Vinci’s designs for mirror-grinding and 
polishing machines. Chapter 2, entitled ‘Single surface working’, describes in helpful 
detail methods for the hand-making of single prisms and lenses. Chapter 3 discusses the 
nature of grinding and polishing and the properties of the ‘polish layer’. 

Chapter 4, entitled ‘ Tools and materials in general use’, is concerned with gauges, 
spherometers, grinding tools, avoidance of contamination by dust during smoothing and 
polishing, use of wax polishers, washing and grading of abrasives, new types of abrasive, 
cements, plasters and varnishes used in optical work. Chapter 5, on ‘ Dioptric substances ’, 
discusses optical glass types, effects of thermal history on refractive index (a subject which 
is taken up again in Chapter 14), the working of artificial crystals and optical plastics. 

Chapters 6 and 7 deal with mass production methods and with the manufacture of 
spectacle lenses, and Chapter 9 with that of prisms and of plane parallel plates. Chapter 8 
{by R. J. Bracey) deals, rather briefly, with the making of microscope lenses. 

Chapter 10 contains an excellent account of the making of non-spherical surfaces, and 
in particular of Schmidt corrector plates. ‘The next chapter, on the testing of optical 
work, is mainly concerned with interference testing, but some account is also given of 
recent work on the Foucault and phase contrast tests (Zernike’s name is misspelt 
throughout Section 227.1 and also in Section 252.3 later on). A description is given of 
the ‘ angle dekkor’ and its use, as well as of some older devices. A full account of the 
Hilger interferometers is reserved for Chapter 12. 

Chapter 13 deals with surface treatments—reflecting and non-reflecting coatings, 
interference filters and the making of graticules ; Chapter 14 describes the testing of 
optical glass and the interesting discovery that, in most practical cases, it is the disturbing 
effect of uneven cooling on the refractive index rather than birefringence troubles which 
necessitate the refinements of modern annealing. 

Chapter 15 discusses the large object-glasses and mirrors needed by astronomers. 
‘The methods of H. Draper, of Sir Howard Grubb and G. W. Ritchey are described, with 
long quotations from the original papers, and the descriptive notes (by Mr. J. V. Thomson) 
on the construction of the 200-inch mirror of the Hale telescope on Mt. Palomar are 
accompanied by many photographic illustrations. Lastly, extracts are given from 
‘Texereau’s important paper on the detection and the effects of small residual imperfections 
in large optical surfaces. The work concludes with a list of reference books and with 
four Appendices. 

The volume is a goldmine of practical information solidly based on experience, while 
its vitality and general level of interest set a new standard for English books on this subject. 
E, H. LINFOOT. 


Contribution a Tétude du coefficient de viscosité des liquides, by J. HUuEeETz. 
Pp. xv+70. (Paris: Publications Scientifiques et Techniques du Ministére 
de l’Air, No. 260, 1952.) 750 fr. 


Contribution a Vétude aérodynamique de l’aile et de Vhélice, by M. MEnarp. 
Pp. xvi+187. (Paris: Publications Scientifiques et ‘Techniques du 
Ministére de |’Air, No. 262, 1952.) 1200 fr. 


Chromophotographie des vibrations dun fluide, by J. M. Bouror. Pp. xviii+77. 
(Paris: Publications Scientifiques et Techniques du Ministére de I’Air, 
No. 264, 1952.) 800 fr. 


Molecular Microwave Spectra Tables, by PAuL KisLuik and Cuartes H. ‘Townes, 
National Bureau of Standards Circular 518. Pp. vit+127. (Washington, 
D.C.: U.S. Department of Commerce, 1952.) 65-c. 
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ABSTRACTS FOR SECTION A 


A Theoretical Attempt to Predict the Excited States of Nuclei in the Neighbourhood 
op - “Fo by M. He L. Pryce, 


' ABSTRACT. *8Pb has a double closed shell of 82 protons and 126 neutrons. The 
classification of the lowest excited states of 2°7Pb, 2°’Tl, 2°Pb, ?°°Bi, which differ by one 
nucleon, should thus follow the single-particle model closely. The attempt is made to 
classify the lowest excited states of 2°®Pb, 2°*T'l, ?°*Hg, 7!°Pb, ?!°Bi and 24°Po, which differ 
by two particles, in terms of the one-particle model, taking the interaction of the two particles 
as a perturbation. Considerable uncertainties in the order of the energy levels remain, but 
some definite conclusions are possible. The existence of an isomeric state of #1°Bi (RaE) 
with very high spin ([=8 ?) is predicted. ‘This fits in with known facts about the alpha- 
particle emitting isomeric state. 

_ The method can also be applied to *°*T', ?°*Bi and the excited states of 2°°Pb which arise 
by lifting one nucleon into the next shell, leaving a hole. It meets with considerable success 
for *°8Tl. Discrepancies with known facts in the case of ?°8Pb suggest a more complex 
picture, according to which the first excited state would have J=0 and even parity. Sucha 
state has never been observed, but its existence is not completely excluded by observation. 


The Angular Distribution of y-Radiation from Aligned Nuclei, by N. R. STEENBERG. 


ABSTRACT. The angular distribution of y-radiation from nuclei aligned at low tempera- 
ture is given as a series of Legendre polynomials, [(@)=1-+L,Ax(T) Px (cos@) the Ax being 
temperature dependent. The formulation is applicable to cascade processes of arbitrary 
multipole orders and nuclear spins. An approximation valid for low degrees of alignment is 
given, as is the angular distribution from completely aligned nuclei. Application is made to 
©°Co. 


A Formulation of Beta-decay Theory for Forbidden Transitions of Arbitrary 
Order: I—WSelection Rules and Energy Spectra, by J. A. SPIERS and 
R. J. BLIN-STOYLE. 


ABSTRACT. The matrix elements of beta-decay theory are expressed in a way which not 
only gives the beta-spectra and angular distributions for arbitrary degree of forbiddenness, 
but also yields very simply the selection rules appropriate to any type of interaction. 
Results for the energy spectra confirm those found by Greuling by inspection ; angular 
distributions and electron—neutrino correlations are considered in a subsequent paper. 


A Formulation of Beta-decay Theory for Forbidden Transitions of Arbitrary 
Order: II—Angular Distributions, by J. A. Sprers and R. J. BLIN-STOYLE. 


ABSTRACT. ‘The method of a previous paper is used to derive general formulae (valid 
for arbitrary Z) for the angular distributions of beta-emissions. 


Studies in Graphite and Related Compounds: I—Electronic Band Structure in 
Graphite, by C. A. CouLson and R. ‘T'ayLor. 


ABSTRACT. ‘The Bloch (tight-binding) approximation is applied to a study of graphite. 
The mean energy per atom and the internuclear distance are calculated in good agreement 
with experiment. Previous accounts are extended by including certain overlap integrals 
and by dealing with a set of parallel layer planes instead of only one plane. ‘The shape 
of the band due to the z-electrons is calculated, and the experimental x-ray emission 
spectrum of graphite is interpreted. It is concluded that the role of o-electrons is greater 
than is often supposed. 
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Studies in Graphite and Related Compounds: I]—Momentum Distribution in 
Graphite, by W. E. DuNcaNsoN and C. A. COULSON. 


ABSTRACT. The space wave functions for the 7-electrons in a graphite layer, as found 
in Part I, are transformed into momentum coordinates, and the radial distribution function, 
including o-electrons and z-electrons is calculated. Using the theory of Jauncey and 
DuMond the shape of the Compton profile in x-ray scattering is predicted. Both this, 
and the momentum distribution, are in excellent agreement with experimentally determined 
values. The introduction of a plausible scale factor in the atomic wave functions would 
make the agreement almost perfect. 


Studies in Graphite and Related Compounds: I11]—-Electronic Band Structure in 
Boron Nitride, by R. TayLor and C. A. COULSON. 


ABSTRACT. An analysis is made of the distribution of z-electrons in the crystal of 
boron nitride. The insulating character of this substance is shown by the fact that these 
electrons completely occupy one band, of width approximately 2-4 ev, separated from the 
next higher band by a gap of about 4-6 ev. The density of states in the band is 
calculated, accurately for a single layer plane, and approximately for the three-dimensional 
crystal. Until the energies of the o-electrons are better known, further progress in 
understanding this substance will be difficult. 


The Diamagnetic Anisotropy of Large Aromatic Systems: IV—The Polyacenes, 
by R. McWEeny. 


ABSTRACT. 'Thez-electron diamagnetic susceptibilities of polyacenes of up to ten rings 
are calculated by a method introduced in an earlier paper: the method, which involves only 
a single numerical integration, is simplified. An alternative method of calculation, which has 
certain attractions and is widely applicable, is sketched in an Appendix. 

The relationship between the diamagnetic anisotropy and the number of rings in the 
molecule is found to be almost linear. 


Strong Coupling in Inelastic Collisions of Electrons with Atoms, by H. S. W. 
Massey and C. B. O. Monr. 


ABSTRACT. ‘The various approximations which have been used to calculate the 
probability of excitation of an atom by impact with an electron whose energy is not greatly 
in excess of the threshold are known to fail under certain conditions. This is especially 
true in the important cases in which the excitation proceeds mainly through electron 
exchange. In these cases the exchange coupling is too strong for approximate methods, 
which assume weak coupling, to be valid. By using a schematic model exact solutions 
for strong exchange coupling are obtained and compared with the results of different 
weak coupling approximations. It is found that the Born—Oppenheimer approximation 
is never reliable at low electron energies but if distortion of the incident and scattered 
waves by the atomic field is allowed for correctly much better results may be obtained. 


Note on the Fluctuation Problem in Cascade Theory, by H. MesseL and 
R. B. Ports. 


ABSTRACT. A general relation between two fundamental functions appearing in 
fluctuation problems in cascade theory is derived. This leads to a direct method of solving 
for the 7th moments of the distribution function. It is also shown that the Janossy 
G-equations are superfluous for the solution of the fluctuation problem. 
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